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Foreword 


It is often said that scientific texts quickly become obsolete. 
Why are the Pauli lectures brought to the public today, when 
some of them were given as long as twenty years ago? The 
reason is simple: Pauli’s way of presenting physics is never 
out of date. His famous article on the foundations of quantum 
mechanics appeared in 1933 in the German encyclopedia 
Handbuch der Physik. Twenty-five years later it reappeared 
practically unchanged in a new edition, whereas most other 
contributions to this encyclopedia had to be completely re- 
written. The reason for this remarkable fact lies in Pauli’s 
style, which is commensurate to the greatness of its subject 
in its clarity and impact. Style in scientific writing is a quality 
that today is on the point of vanishing. The pressure of fast 
publication is so great that people rush into print with hur- 
riedly written papers and books that show little concern for 
careful formulation of ideas. Mathematical and instrumental 
techniques have become complicated and difficult; today most 
of the effort of writing and learning is devoted to the acquisi- 
tion of these techniques instead of insight into important 
concepts. Essential ideas of physics are often lost in the dense 
forest of mathematical reasoning. This situation need not be 
so. Pauli’s lectures show how physical ideas can be presented 
clearly and in good mathematical form, without being hidden 
in formalistic expertise. 

Pauli was not an accomplished lecturer in the technical sense 
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of the word. It was often difficult to follow his courses. But 
when the sequence of his thoughts and the structure of his 
logic become apparent, the attentive follower is left with a new 
and deeper knowledge of essential concepts and with a clearer 
insight into the splendid architecture of reason, which is theo- 
retical physics. The value of the lecture notes is not diminished 
by the fact that they were written not by him but by some of 
his collaborators. They bear the mark of the master in their 
conceptual structure and their mathematical rigidity. Only 
here and there does one miss words and comments of the mas- 
ter. Neither does one notice the passing of time in his lectures, 
with the sole exception of the lectures on field quantization, 
in which some concepts are formulated in a way that may 
appear old-fashioned to some today. But even these lectures 
should be of use to modern students because of their compact- 
ness and their direct approach to the central problems. 

May this volume serve as an example of how the concepts of 
theoretical physics were conceived and taught by one of the 
great men who created them. 


Victor F. Weisskopf 


Cambridge, Massachusetts 


Preface 


“Statistische Mechanik” was the first of the six volumes of 
lecture notes of Pauli published by the “Verein der Mathe- 
matiker und Physiker an der E.T.H. Ziirich.” After many years 
of lecturing, Pauli’s personal notes had started to disintegrate. 
So he asked his assistant, who at that time was the late M. R. 
Schafroth, to take notes for his personal use. The result was a 
short compendium with almost no text. But the students asked 
Pauli for the permission to publish the notes for their use, 
which they did in 1947. It was indeed far from easy for an un- 
initiated student to follow a course given by Pauli. 

But Schafroth’s notes did not only reflect the original pur- 
pose of a compendium for Pauli’s personal use. They were also 
typical of Schafroth’s style of talking physics. It is the same 
style as that of the Selected Topics in Field Quantization, in 
this series, which was also based on notes by Schafroth. 

These particular circumstances of the origin of the present 
lectures of course posed many problems both to the translators 
and to the editor. Phrases had to be inserted and some deriva- 
tions needed changes and amplifications. The result, I believe, 
is a concise course on statistical mechanics giving today’s stu- 
dent a quick, if not easy, introduction to the subject, centered 
on the historic development of the basic ideas and on the logi- 
cal structure of the theory. 

Since these lectures were given by Pauli important advances 


have been made in the rigorous solution of simple models of 
ix 
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statistical mechanics which have given detailed insight into 
the phenomenon of phase transitions. The first and most spec- 
tacular achievement in this field, Onsager’s solution of the 
two-dimensional Ising model, although published three years 
before the German original of these lectures, had not been dis- 
cussed by Pauli. Most likely at that time Pauli considered this 
problem too advanced and technical a subject for the level of 
this course. And, in fact, since then the technique of deriving 
Onsager’s result has been considerably simplified. A good in- 
troduction to this field and, more generally, to the problem of 
phase transitions can be found in the new book by H. E. Stan- 
ley, Phase Transitions and Critical Phenomena (Oxford Uni- 
versity Press, New York, 1971). 


Charles P. Enz 


Geneva, 19 November 1971 


Pauli Lectures on Physics: 
Statistical Mechanics 


Chapter 1. Stosszahlansatz } 


1. CONCEPTS FROM THE ELEMENTARY KINETIC THEORY OF GASES 


Let f(v)d*v be the number of molecules with velocity 
contained in d?v, and let n= | f(v)d?v be the number of 
molecules per cm*. The function f(v) can also depend on x. 
Then, we have 


| [tee x)d3v d3x =| m(xare =—N, 


where WV is the total number of molecules. 
The pressure tensor is defined by 


Dik = M]0;%,f dv, Pik = Pri = NMViU; , 


where m is the mass of the molecules. The mean velocity 
is defined by 


If there is molecular flow, c 40, we define 


u=—v—cC. 


1 The ‘“‘Stosszahlansatz’”’ is an assumption related to the calculation of the 
number of collisions between the molecules of a gas. 
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Then, strictly speaking, the pressure tensor is defined as 
if — NMU; Uy, = mM U,;U,f(u + c)dsu. 


Since u=0, we have 


UV, V_ = C, Cy Uz Uy « 


Example: Consider a homogeneous, isotropic velocity dis- 
tribution: 
c,=0, f(v) = F(r) , v=|v|, 
Pr=P0n, pP=nmv?i=nmei=nm?. 
Thus, 
p= nmr. 


The ideal gas laws follow from this equation: Since 
inmv = U =average kinetic energy of the molecules 
per cm’, the equation can be rewritten as 


p=8U. 


2. COLLISION LAWS 


We consider collisions between two masses m and M 
which move with velocities v and V, respectively. We de- 
fine the relative velocity by 


w=v—V 


and the velocity of the center of mass by 


pet mv+ MV 
~ m+M 
Solving for » and V, we obtain 
M 
re Ot ae Mim” 
[2.1] 
|| ee 
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These equations, combined with the definition of the kinetic 
energy, E,,,= $mv?+ iMV?, give 


M 
m + eee mM. [2.2] 


Te 2 
" 2 2M+m” 


If primes are used to denote quantities which refer to times 
after the collision, then the conservation laws of momen- 
tum and energy may be written as 


P= P>U'=U 


; [2.3] 
Eyn = Exin 
Substituting U’=U into E,,,=F,,, gives 
ipa ote aw) [2.4] 


The conservation laws do not suffice for determining the 
variables after the collision. There are two pieces of infor- 
mation still lacking, which we can determine by means of 
special models. 


a. Billiard ball model 


We idealize the molecules by considering them to be 
rigid, elastic spheres. The collision is specified by a unit 
vector n which points from the center of m toward the 


Figure 2.1 
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center of M at the time of collision. We then have 
f / 
10, = “18 tO — oe [2.5] 


In place of nm we can introduce two angles: 6= angle 
between w and n, and e=angle between the plane deter- 
mined by w and n and a fixed plane containing w. Then, 
the scattering angle is given by y=2— 26. If we vary 
the initial conditions a little, so that U is contained in d°U 
and w is contained in d3w, then U’ will be contained in 
d?U' and w’ will be contained in d*w’. The differentials 
are related by 


d?3U0=d3U' (from Eq. [2.3]) 
d3w = dw’ (from Eq. [2.5]). 


Therefore, 
awd? = dtw'dtU’ . [2.6] 


Furthermore, because of Eq. [2.1], we have 


ail 2. O(U;, w;) 
OV,v) avn 


Using Eqs. [2.6] and [2.7], we obtain 
~ a8V dev = a3V' deo’. [2.8] 


=—!xixl=— [2.7] 


If we allow nto vary within the solid angle d?A= sin@ d0 de, 
then, because there is no d?A’, we have 


d?V d3v d2A = d?V' d8v'd2/ . [2.9] 
This is a collision invariant which we shall find useful later. 


b. General central force model 


We assume that the molecules are mass points which 
repel each other according to some central force law. We 
can draw a diagram analogous to that of Fig. 2.1, in which 
the straight lines along w and w’, instead of being the 
trajectory of the center of mass, are the asymptotes to 
the trajectory. 
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ia) 


Figure 2,2 


Then we can introduce the unit vector n which is par- 
allel to the bisector of the angle between w and w’ in the 
plane of the trajectory, and which points in the direction 
from m toward M. We can also introduce the impact pa- 
rameter b, the scattering angle y= y(|w|, b), and the angle 
of the plane of the trajectory «. As in Fig. 2.1, we define o 
by o=b/sin6, where 6=(x—vy)/2. Then Eq. [2.5] still 
holds; therefore, Eqs. [2.6], [2.8], and [2.9] are also valid, 
and d*vd°V is a collision invariant also in this case. 

We introduce two additional concepts: the inverse col- 
lision and the opposite collision. They are characterized 
by the parameters in the table below. 


Original Inverse Opposite 
Collision Collision Collision 
v ,V;0,¢ —v',—V'; Oe ara initia 
won —U',—w';n U',w'; —n 
ES —v,—-V v iV | snat 
U'",w' —U,-—w U,w 


3. CHANGE OF THE DISTRIBUTION FUNCTION DUE TO COLLISIONS 


We assume that there are two types of molecules: type 1 
is characterized by f and v, and type 2 is characterized 
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Figure 3.1 


by Fand V. A decrease in the number of molecules of type 1 
in d*v by collisions arises from a) collisions with molecules 
of type 2, and b) collisions with molecules of type 1. The 
number of collisions per second, characterized by b and e, 
between molecules of type 2 whose velocities are in d*V 
and molecules of type 1 whose velocities are in d*y is 


wh db def(v) F(V) dsvd8V . 
(see Fig. 3.1). This implies the following (negative) rate 
of change of the number of molecules of type 1 in d*v: 


aot = =| [farvaracwose)rr)arv. 
The rate of change due to the decrease of the number of 
molecules in d°v because of collisions between like mole- 


cules (type 1) is 


dev t= —[f farvasdewd jo) f)are. 
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Similar consideration of the increase of the number of mol- 
ecules in d°v due to collisions leads to 


av = +|{faveas dewb f(v’) F(V’) d3v' 
+{ffarvras dewb f(v') f(V') dv’. 
In this expression, the given quantities are » and V rather 
than v’ and V’. It would be more precise to write 
Do —2 (v,.V 0,2) and y= (v, 1b, 6) « 
By application of Eq. [2.8] it follows that 


— f= [[fmraraearr (poy m( — fv!) FV’) 


+] {[uvaracav toi) — 104). 
The impact parameter is a function of 6 and w, b=5(6, w). 
Therefore, 


bdbde = b (5). d6de= Qda?A, 


ob) 1 
ae elses : 


For rigid spheres, @=o*?cos§. (In what follows, we use Q 
for M-M collisions, Q for M-m collisions, and q for m-m 
collisions.) 

Finally, we have 


meee) =| fn@araarv [joy r — f(v') F(V'))] 
£ i wq@ Aa fv) (V) — fv (V 1, 


where 


and 
EAP) [ fudaraarot po 7h) — fo) 
3 i wQ a2Ad*e[F(v) FV) — F(v') F(V')) . 
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Remark: In deriving this equation, we have introduced 
a fundamental hypothesis which is in addition to the as- 
sumption of central forces only; this hypothesis is called 
the ‘‘Stosszahlansatz.’? We have assumed that the density 
of molecules in the cylindrical volume element bdbde is 
the same as it is in the remaining gas. Since, however, 
there will be density fluctuations, of/ot and OF /ot will also 
fluctuate. By neglecting these fluctuations, we actually 
have not calculated of/dt; instead, we have caleulated Af/At. 
This one-sided average also means that our formula picks 
out a specific direction in time; this is in contrast to the 
laws which were used in the derivation. 


4. STATIONARY DISTRIBUTIONS 
Case 1. Distribution functions independent of position 


We define the Boltzmann H-function, 
i, =| ftogsare + |FlogFd?V, [4.1] 


and a related function, 


H = [fog f—1) a 


+] P(og¥—1)a°7 = H,— Nm— Ny » [4.2] 
Because of number conservation, it follows that 
aH _ dH, 
dt dt 
=ftogs fare + log F oF av 


—[ | fu@arzaroarv fog f(v) + log F(V)] 
x[f(v) FV) —f(v) FW) 
x { | | wq d?Ad*” dV [log f(v)][f(v) (VF) — Hv) (0) 


-{ ill wQ dA d*v d°V [log F(V)][F(v ene) 
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We have 


[ [foca:zare ar tog (oye) 10) 


1 

= 5] | froraez ar ary [og fo) + tog 1710) HP) 
1 

=} | | i oq 2A. dv dV flog [f(v) {(V)} f(v) #(V) 


(b) Because the molecules are identical, interchanging v 
and V also interchanges v’ and V’: 


| fom d?/ d*v d°V [log f(v)] f(w') (VV 


=51| [og d*4 d°v d°V flog [f(v) ((V)} flo’) (P). 


(c) With these results, we obtain for the second integral 
in the expression for dH/dt, for ial 


—5 | [[ocaravarr foster} (ie (v') (V)}. 


(d) Because of the properties of opposite collisions, this 
integral can be written as 


=| | {roaaraare a7 flog f(a) (¥)]— lost fo ie 
x {f(v)f(V) — fv) f(V')} . 
Finally, it follows that 


= = 3 | { fo@araare d°V {log[f(v) F(V)] — log[f(v’) FV’) }} 
xo) FV) — fv FV) 
ey 7] feoe d2/A d3v d*V {log[f(v) f(V)] — log[f(v’) f(V’’) }} 
x {f(v) f(V) — flv’) fV} 
-s a) I | wQ a?/. dev dV flog [F(v) F(V)]— log[F(v') F(V')}} 
x {F(v) FV) — F(v') F(V)} . 
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Since (we—y)(loga—logy)>0, for gw, 
=0, for 7=y, 


we therefore obtain the Boltzmann H-theorem: 


di <0 [4.3] 
We are seeking stationary conditions, that is, distribu- 
tions f and F for which of/ot=0 and oF /ot=90. From this 
it follows that dH/dt=0, which is therefore a necessary 
condition for a stationary distribution. Because dH /dt=—0, 
we have 


fv) (V) = f(r’) {V') [4.4] 
F(v) F(V) = F(v') F(V’) [4.5] 
f(v) FV) = f(v') FV’). [4.6] 


These relations must hold for all vw’ and V’ which are al- 
lowed for the values of v and V; that is, they must hold 
for all vw’ and V’ which satisfy the energy and momentum 
conservation laws: 


p= mv+ mV = mv'+ mV'= p', [4.7a] 

P= Mv+ MV= Mv'+ MV'=FP’, [4.8a] 

P= mv+ MV= mv'+ MV'=P', [4.9a] 
m m Mie, M_, 

e=yet sage ts =e, [4.7b] 
M M M sae : 

BS + > ast ee ee [4.8b] 

Rigas Ops Mt Pes Le UM iS oR 

E 5 0) erred oo ae gee ee [4.9b] 


However, from this it follows that of/ot=0 and oF /ot=0; 
that is, dH/dt=0 is also a sufficient condition for a sta- 
tionary distribution. If we substitute 


logf=gp and logF=@, 
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then Eqs. [4.4]-[4.6] become 


pv) + OV) = p(v')+ oh), [4.4a] 

P(v) + OV) = Div') + DV’), [4.5a] 
and 

pv) + BV) = ov’) + OV’), [4.6a] 


with Eqs. [4.7a]-[4.9b] as auxiliary conditions. 
If we treat the problem with Lagrange multipliers, we 
obtain 


E ia eee @m| ‘dv +[...]-dV=[...]-dv’+[...]-dV". [4.4b] 


The quantity in each bracket must equal zero: 


oy, 
ae a'mv + Bm =0. 
The solution is 

f = A’ exp[— «dmv? —B-mv]. 


If we define c=—B/a’ and a'm/2=a, A’ exp(ac?)=A 
then 


f= Aexp[— a(v—c)?]. 


This is a Maxwell distribution with a superimposed constant 
molecular flow c. 
Similarly, from Eq. [4.5a] we obtain 


F= A'exp[—a’4 MV?—b-MPV). 


Since, from Eq. [4.6], fF can only depend on £ and P, 
it follows that a’=a«’ and b=8, which means that c is 
the same in both cases. Therefore, the two gases may not 
move relative to one another. Thus, it follows that the 
most general stationary distribution independent of posi- 
tion is the Maxwell distribution: 


f = Aexp[— a(v—e)?], 
F = Aexp[— a(v—c)?]. 
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Case 2. Dependence of f on DOSES and effects of external 
forces 


These two factors give rise to additional terms in the 
expression for of/dt. 


(a) Dependence on position. From 
f(t+ dt, x, v) = f(t, x— vd, v) 


follows 


a; aa. 


oy Ax vat. 


Thus, the contribution to the time derivative of f due to 
motion of the molecules is 


(b) Effect of external forces. In this case 


hence 


f(t+ dt, x, r)=t(t x, v— Tat) : 
m 


Thus, the contribution to the time derivative is 


With these results we obtain the complete basic equation 
of the kinetic theory of a one-component gas [A-1] ?: 


0 of K 
at Se Ox vt =| favaratte’ IV) — flv) fV)]wq . 


* Comments [A-1]-[A-3] appear in the Appendix on pp. 111-113. 
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The generalization to two- and more-component gases is 
immediate. 

In order to find stationary distributions, let us assume 
that K is a conservative field: 


OL pot 


a Ox - 


From the assumption that the distribution is stationary, 
it follows that of/Ot=0 and that the change of f by col- 
lisions must vanish; that is, the right-hand side must 
vanish. Therefore, 

(1) f is a Maxwell distribution with position-dependent 
constants, and 


(2) mv 01 — OB oF _ og. 


This is satisfied if f depends only on F,,, and £,,,: 


of of OF 50 
Ox OF Ox ’ 
of of 
dv on? 
hence 
of OBye Of 


oo Ox ov 


Thus, we obtain a stationary distribution if we replace 
L,.= ¢mv? by H= EL, +£,,, in the Maxwell distribution. 
The resulting distribution is called the Maxwell-Boltzmann 
distribution: 


{(v, x) —"A! exp |- a’ ¢ v? + Ba())| ; 


Remark concerning the exact derivation of this formula. We 
must begin with a generalized H-theorem in which the 
position-dependent H-function is integrated over the gas 
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volume @: 


KH = aaa j= fara (flog #—f), 


dx? Of 
3 3 Foy ul 
ae a of vlogi a, 
4 
sig lool log 7 SE v = ppm Fe K 
ov m 


= 5 —ate{ aston — 1) — [Br eat f). 


Here, J<0 and 2 is the boundary of G. The integral over 
the infinite sphere 2 in v-space clearly vanishes. 

The second integral can be made to vanish if we assume 
that consists of ideally reflecting walls, or that the gas 
is not bounded, so that X can be placed at infinity. Then, 


dH? ceils ‘ 
a J =fa ofa ofa vfa A 
x [fw ) (VV) — fv) f(V) qu log f <0. 


For stationary distributions we have of/ot = 0, which implies 
d#/dt=0. Again, this is only possible for 


Fv')f(V') = fv) f(V) 
with specified auxiliary conditions. The result is 
{= Aexp[— a(v—c)*], 


where A, a, and ¢ are functions of x. This satisfies the 
basic equation and the condition for a stationary distri- 
bution, df/ot=90, only if [A-1] 


of K of dx of dw 
oes m Ox ae “onee at a 


af 
Ox 
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that is, only if f is a time-independent integral of the 
equations of motion of the molecule. From this, we can ob- 
tain the following most general distribution function [A-1]: 


Bien Byes) 20-e)] 


rafal 


m 
Here, (1) «constant and A=constant, 
(2) K must have a potential Z,,,, and 


(3) c= c®+wx x corresponds to the motion of a rigid 
body, so that the velocity vectors are everywhere 
tangential to the equipotential surfaces. 


5. TRANSPORT PHENOMENA 


The fundamental equation also provides the exact basis 
for treating transport phenomena. 


a. Viscosity 


We neglect external forces: 


af of =| farvaraongigo (0) — flv) {(V)). 


The pressure tensor, 
Dik =| M0; Vz f dv , 
is to be calculated. As an example, we consider the case 


of a gas between two plane plates, one stationary and one 
moving. As a zeroth-order approximation, we can set 


¢,(z) = C— (from the macroscopic theory); 


then, 
fo(v) =A exp[— a(v — ¢)?] 


ayers -« [(o-- of) +e +e]| 
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However, this is not a solution of the fundamental equa- 
tion, because of the term (éf/éz)v, on the left-hand side. 


Figure 5.1 


We obtain a further approximation by writing 


f=f,[1+ uu, B(ut) +...], 
with 


&, =,; 4, =O): 


Ral a 
— 


The fundamental equation yields 


S- | JaeA aera tf rfol UY [ue we Bow) 
+ Cel BIT") — fait Vi eaeBuw) + ER ep . 


It is easily verified that the dependence on the direction 
of the velocity which results from the spherical symmetry 
of the forces is identically satistied. What remains is a 
complicated integral equation for B(w?). (See S. CHAPMAN 
and T. G. CowWLING, Mathematical Theory of Non-Uniform 
Gases.) The shearing stress can be calculated in terms 
of B(u?): 


V2= m| ut o(s) B(u?) ase. 
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Without calculation, it is clear that the result must be 


C OC, 
Bie hee al oe ees 


Remark: Our approximation makes the shearing stress 
linear. By using the concept of the mean free path 1, which, 
of course, is implicitly contained in the integral equation, 
we can, on the basis of the elementary theory, suppose that 
the approximation will be good only if J is small compared 
to the dimensions of the vessel. The higher-order terms, 
which must be considered if this condition is not satisfied, 
result in deviations from the macroscopic theory. 


b. Heat conduction 


The heat conduction problem can be solved with the 
more general expression [A-1] 


f(t) = fo() (1+ 5 (ian — 4dut) £2 Beas) 


02, 
oT 
4. >t ap 4. 4 ’ 


where JT = m/(2ka) is the temperature and k Boltzmann’s 
constant. 

We now compare the results of the elementary theory 
with those of the rigorous theory of Chapman. With the 
elementary theory, the viscosity coefficient is given by 


where y~1 is an undetermined factor; the heat conduc- 
tion coefficient is given by x=6yc,, where 6 is another 
factor of order one. With the rigorous theory, we obtain 


which means 
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we also obtain 
% = 2.522 Xe, , 
which implies 
& = Baz . 


Note: (1) The form of the law of force between the mol- 
ecules does not affect the lack of dependence of 7 on pres- 
sure; it does, however, affect the temperature dependence 
of 7. 

(2) The fact that the only tensor which enters is one 
whose trace is zero, i.e., uu.—46,u?, Is a result of the 
spherical symmetry of the forces, and it leads to the con- 
sequence that the pressure tensor is zero for spherically 
symmetric velocity distributions. In other words, the two 
constants, which are possible in principle because we cal- 
culate with tensors, reduce to one constant: 


Oc; , OC 
pa=—n| (ae +5¢) iba =e sch Lia) 


6. MEANING OF THE H-THEOREM. TEMPERATURE 


From Section 4 we have 
te [evarv(fiog/—f = #,—N. 
We now specialize to a Maxwell distribution: 


3 
f=Aexp[— av?], [raro—a (2) =n, 
aw 
so that 


Then, 
= [arearop tog 4 — av?—1) = N (log A — av?— 1) 


= N (logn + 3 log a — $loga — $) 
= N (log N— log V+ $loga— %logm— 8). 
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Furthermore, 


2 
waa N (- log V— 3 log H+ [log +8 og AS) oe | 
a a 
= N(— log V— 3log #+ constant) . 


For the thermodynamic entropy S we have 
1. S is only determined to within an additive constant, 
2. p=T(oS/oV),, 1/T = (0S/e#),. 
For one mole of an ideal gas, 


S= C,logH+ RlogV-+ constant. 


On the other hand, for one mole of a monatomic gas, 
3L 
KH = — “9 log E— Dlog V+ constant , 


where Z. = Loschmidt’s number (Avogadro’s number). Com- 
paring, with R/L=k, we obtain 
S=— k#-+ constant, 


and O,=3, which confirms a previous result (see end 
of Sec. 1). From this it follows that 


es-.(-. 
es V 


oH vu oH 2H m’ 
where 
m m 
Sa? tae 


Remark: In thermodynamics, temperature is the primary 
quantity; from this, with the help of general assumptions 
(the first two laws), the concept of entropy is derived. 
On the other hand, in statistical mechanics a function # 
can be defined which has the properties of entropy ; from # a 
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temperature can then be defined as a secondary quantity, 


ox\ 1 

On}, kr’ 
as was done above for a monatomic ideal gas. Using the 
concept of temperature, # can be written as 


Qak 5 
—-H= W (log 5 + Hog T+ §log==" +5) 


We must not ascribe too much meaning, in the sense of 
Nernst’s theorem, to the entropy constants which appear 
here, because ¥ is not invariantly defined with respect to 
dimensions. Indeed, fd*vd?a# is dimensionless, but f is not; 
therefore, { {flog fd®vd%« is only defined to within an addi- 
tive dimension-dependent constant, and the indeterminacy 
of the entropy constants remains. This question is only 
resolved in quantum statistics. 


7. STATISTICS OF AN IDEAL GAS 


We divide (v, x)-space into cells, each of which is of 
constant volume w, and which we number with an index k; 


N,, = f d?v d3zx 


is the time-varying occupation number of the kth cell. 
Since the N;, change discontinuously, we only consider finite 
time intervals rt. 

We call the set of all M, a macrostate of the gas. Of 
course, we always have 


SM=N, and > N,e, = HE = constant, 
k k 


with e,=mv;/2 for our case of a monatomic ideal gas. 
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A microstate of the gas is defined to be a set of numbers 
which specify in which cell each atom is located: 


s =number labeling the atom, 
k, =index of the cell in which atom s is located, 


§ 


(k,) = microstate. 
We have 


Se 


s=1 


The macrostate is uniquely determined by the microstate, 
although the converse is not true. For every macrostate 
there are very many microstates, i.e., as many as the num- 
ber of ways in which N elements can be ordered into 
groups of N,, N2, ... identical elements. This ‘‘multinomial 
coefficient”’ is well known: 


M! MN! 
ANN. Lee 
k 


Boltzmann’s Fundamental Hypothesis: All microstates are 
equally probable. 

Discussion: We must consider “probability”? as relative 
(a posteriori) frequency of occurrence, but in what kind 
of ensemble? The possibilities are 
1. Time ensemble of a gas, consisting of discrete times 
out of the time development of the gas, 

2. Statistical ensemble, consisting of very many gases at 
one instant in time. 

Boltzmann assumed that these two ensembles give the 
same results. The time ensemble is determined mechani- 
cally, and it can be expected that this hypothesis will not 
be satisfied for certain special initial states. The weight 
of these exceptional initial states should, however, be van- 
ishingly small compared with that of the other initial states. 

For the validity of the Boltzmann hypothesis, it is nec- 
essary that there be no integral of the equations of motion 
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other than the energy; this is immediately clear. Sufficient 
conditions are less trivial. 

The Boltzmann hypothesis is closely related to Gibbs's 
ergodic hypothesis, which apparently cannot be proved 
within the framework of classical mechanies. 

If we accept the Boltzmann hypothesis, then it follows 
that the relative probability of the macrostate (¥,) is 

Ni! 
W —— 


~ [ae 
k 
For ¥,>>1, using Stirling’s formula we obtain 
log W = log N!— > NW, (log ¥,—1). 
k 
a. Relation to # 


H -| it (log f — 1) d*v dSxv 


N, 
= N, ] = 
54: (lor —3) 
= > Ne (log Ne — 1)— N log ; 


hence, 
log W = log N!— #— N logw. 


The first and third terms are independent of the distribu- 
tion; they compensate for the different dimension of 


b. Most probable distribution 


Using Lagrange multipliers, we obtain the following result 
for the extremum: 


Clog W . 
a: —a&:+B=0, N,=Aexp[—ae]. 


That is, the Maxwell-Boltzmann energy distribution is the 
most frequent macrostate. 
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c. Theory of fluctuations in ideal gases 
We consider states in the neighborhood of the Maxwell 
distribution: 
Nz = A exp[— ae; ], 
M, = N2+4,, 2 4:=9, x &%4r=0, 


olog W 
OM: ) 4. 


log W = log W, + oH 
k 


o? log W 
aS 4 2 (Syeane), 4940+ sees 


The linear term vanishes because of the auxiliary conditions, 
and the quadratic term becomes 


onog Wl 3 
cagoN, mm *** 
Thus, 
ta iy 
log W = log W,— 4 > — 4i, W=W, [| exp = =weroll* 
~ Ni k 2N;, 


We can draw the following conclusions: 

1. The Maxwell distribution actually represents a maxi- 
mum of W. 

2. For density fluctuations, we have 


“A2 70 
pV ae 


d. Discussion of the H-theorem 


After # has reached its minimum, it will not keep this 
value indefinitely ; rather, fluctuations will occur which are 
more infrequent the larger they are. If a large fluctuation 
occurs at time t, it will very seldom happen that an even 
larger one will occur at time t+, or that a larger one has 
occurred at time t—vr. Because of this, time-symmetry is 
reestablished. What we have calculated by means of the 
Stosszahlansatz is not valid for a single gas; it is valid 
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only for a statistical ensemble. With a statistical ensemble, 
the values of #, necessarily discrete for a single gas, are 
replaced by a continuous distribution of values which sat- 
isfies the H-theorem. However, when we thermodynami- 
cally observe irreversibility, the situation is always such 
that the distribution of positions and velocities approaches 
the equilibrium distribution with overwhelming probability, 
although small fluctuations continually occur. 


e. Generalization of the concept of entropy 


The statistical view also permits us to formulate a defi- 
nition of entropy for nonequilibrium states. For two states, 
1 and 2, we have 


leaving the additive constant unspecified, we obtain 
S=klogw. 


Because of the logarithm, and because the probabilities of 
independent states multiply, the additivity of entropy is 
maintained. ~- 


Chapter 2. General Statistical Mechanics 


In the following, we no longer restrict ourselves to ideal 
gases. We seek a general scheme with which to derive the 
thermodynamic functions of state of a mechanical system 
from its mechanical structure. 


8. PHASE SPACE AND LIOUVILLE’S THEOREM 


We consider a mechanical system with NV degrees of free- 
dom which is described by its equations of motion in ca- 
nonical form: 


With specified initial values g? and p?, the integration of 
these equations gives the time development of the system, 
geometrically represented by a curve in the 2N-dimensional 
phase space, the orbit in phase space: 

qi (Gey Pe» 0) = qe 


0 


Di (Gey De, 0) =D}. 


Qs = 4 (Gey Dey t) 
Pi=D: (Gey Pes t) 


with | 


Further, from the equations of motion it follows that 


yt p oH oH 
Pr dp, oat of 


HSE _S td (tes 


25 


26 GENERAL STATISTICAL MECHANICS | Chap. 2 


In case H does not depend on time explicitly, then 0H /ot=0, 
and 
H (p,q) = EH = constant. 

This means that for a time-independent H the orbit in 
phase space lies completely on an energy hypersurface. 
Because the orbits in phase space are uniquely determined 
by the differential equations and the initial conditions, 
two orbits can never intersect. 


Liouville’s theorem 


The volume of a volume element in phase space does not 
change in the course of time if each of its points traces out 
the orbit in phase space determined by the equations of 
motion; that is, 

O(p, 9) 
Ht) — sea 
We can prove this in two steps: 
1. First we show that 


[oD(t, to) /Ot],_,, = 0 


Defining «=t—t, we can write 


ap, 34, 
=1te 3 (4 a +... 


t i}: t=t, 
Hence, 
OD(t, =] _ & rE 
| ot t=t 7 2 Opt r = =t, 
oH oH 
“S(sarent ag? 7. 


oH oH 
= (-shraa* auton) ~° 
2. From the multiplication law for Jacobians, we have 
OD(t, to) _ OD(t, tr) 
aa ot D(t, , to) a 
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Letting t, approach t, we then obtain 


OD(t, te) [OD(t, t,) 
ot =| ot |_Pim). 
Hence 
OD (t, to) 
ot 


= 0. | Dt 7,)— constant . 


Since D(t,, t.)=1, we obtain the desired result, D(t, t))=1. 


9. MICROCANONICAL ENSEMBLE 


We consider statistical ensembles with density o(p, q; t) 
in phase space; that is, in the volume element d*”Q there. 
are o(p,q;t)d’*Q ensemble points. As in hydrodynamics, 
let o/ot represent differentiation at a fixed point, and 
let D/Dt represent differentiation along the mechanical 
orbit in phase space: 


Deas aC gn 
pina 3 (bag tap) 
By definition, 
18) 
pi ee" 2) = 9. 


Furthermore, according to Liouville’s theorem, 


D 
SON) . 
Dit 


Therefore, 


o 


na, 


t 


which may be written as 


Og » 0 1s st) = 
f+ 3 (use +ieae =0, 
Oe OH 0g 0 =<) 
ve eee ee () 
au 2 ES Ogr Ox Ox 
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In terms of Poisson brackets, the statement is 
0@ 
+[H,e]=0. 


If the distribution is stationary, we must have d0e/ot= 0, or 
[H,e]=0. 


This is exactly the condition for a time-independent integral 
of the equations of motion. Thus, a stationary density 
distribution may only depend on time-independent integrals 
of the equations of motion. 

The simplest case of a stationary distribution is when @ 
depends only on the energy: e=o0(H). This is trivially a 
stationary distribution and, as is certainly not easy to 
show, it is also the only such distribution which is regular 
(i.e., the only one whose dependence on p and gq is not too 
discontinuous). Therefore, on purely physical grounds one 
can expect this case to include all reasonable densities. 

Especially interesting are those densities which only de- 
pend on energy and which describe closed systems: 


constant, HK<H<H-+dE#, 


— HE eS 
ia 0, otherwise. 


Such a distribution is called a microcanonical distribution. 
For the systems considered here, the energy shell is gen- 
erally a closed, or at least finite, hypersurface. For this 
reason, the volume of the energy shell is finite: 


Jaro = (E)ak. 


E<H<E4E 
In this connection, there are essentially two points of 

view—that of Boltzmann and that of Gibbs. 

a. Boltzmann’s point of view 


In principle it would suffice to consider one system and 
not a statistical ensemble. What is observed macroscopi- 
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cally are time averages of the form 


= lim ni | q) dt. 


This view is based on the assumption that, apart from a 
vanishingly small number of exceptions, the initial con- 
ditions of a mechanical system do not influence these aver- 
ages. This point of view is physically very satisfying, except 
that such time averages can never be computed [A-2]. For 
this reason the ergodic hypothesis is introduced. The ergodic 
hypothesis states that time averages are identical with sta- 
tistical averages over a microcanonical ensemble for reason- 
able functions f, except for a number of initial conditions 
whose importance compared with that of all other initial 
conditions is vanishingly small. 

The ergodic hypothesis can be viewed as a generalization 
of the fundamental hypothesis of Boltzmann (Sec. 7). His- 
torically, the ergodic hypothesis was the assumption that the 
orbit in phase space comes arbitrarily close to every point 
in phase space during the course of time. Our assumption 
is stronger. As an example, consider the case in which 
f=1 in a small region G, and f=0 outside of G. Then, 
we find that 


time average = fraction of time the system spends in G, 


and 
ensemble average = volume of G. 


This means that the fraction of time the system spends 
in @ is proportional to the volume of G. The original 
ergodic hypothesis does not go as far as this statement. 
The ergodic hypothesis, which probably cannot be proved 
classically, can be shown to be satisfied with overwhelming 
probability in quantum theory. 
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b. Gibbs’s point of view 


Neither the time ensemble nor the completeness of the 
theory is essential. Gibbs simply assumes the canonical 
ensemble, which yields physically satisfactory results. In 
spite of the difficulties which are connected with the point 
of view of Boltzmann, Gibbs’s concept is physically much 
less satisfying. For that reason, we shall adopt the point 
of view of Boltzmann. 

The average value of an arbitrary quantity f(p, q) over 
a microcanonical ensemble is 


te wie? 


ao 


In this formula, m means integration over the energy shell 
(see p. 28) 


E<H<E+dk. 


The integrals can be transformed by means of the relation 
fa27¥ Q= Fp fa27Q 
dH ; 
E<H<F+dF QUE) 
Here 2(E£) is the volume 
Q(E) = | a7Q., 
H<E 


Let «; be one of the p, or q,; then we have 


= JeteH| /Ow,) @2*Q daH(d/dE) Bk (OH /Ox,)d2¥ Q 


On, fa?* 2 7 o(H)d# 
d/dE) J ax,(OH /Ox,) d2* Q 
QUE) 
w(£) 
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Since 0H/Ox,=0, we have 


fas _ |e ee 
ard = tai | ay : 


Q(B) QB) 


Performing a partial integration, we obtain 


— [ag @—Hyara=+[n— mara, 


QUE) Q(B) 
Then, 
af Hara 4. fon—mara 
Q(B) a 
=|d%¥O04+(E—H)| =Q(B). 
| H=E 

Q(B) 

Finally, 

0H QE) QE) 1 


"3, w(E) Q(B) (d/dB)log Q(B)’ 


Therefore, this average value is independent of 7. 
quantity 


2 =lag 2 = 2 
k 
has the properties of entropy (see Part ¢, p. 33). 
particular, 
az) 11 
dE}, O kT 
Therefore, 
oH 
wv; én, =60. 


Corollary 1: Let 2;=p,; then 


f= 
DP: op, = D):4:- 


31 


The 


In 
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If, as is usually the case, H is of the form 

H = Eyo(q) + Lun(P,g), where Fyn= 2, %e(Q) PsP , 
then 


Ey Oe OH 
os _ 3 oe 
Op; 2 Par ms 


> Pi 
This allows us to define 
$p,0H/Op; 


as the kinetic energy of the ith degree of freedom. From 
this, the equipartition theorem follows: 
The average kinetic energy per degree of freedom is $kT. 


Corollary 2: Let #;=@q;; then 


This is the virial theorem [A-3]: 


—q;,K,=kT. 


In addition, we can derive some support for the ergodic 
hypothesis. It can be trivially shown that the following 
time average vanishes: 

t 


d 
dt (PiQi) = 0. 


And, from the above, the statistical average over a micro- 
canonical ensemble is 


ai Pidi) = UP Pi = q Ba, Pig = 


If this were not so, then the ergodic hypothesis would be 
severely shaken. 
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Remark 1: If, in analogy with the above, we calculate 


oH 
Ce? 


for 1k, then the partial integration contributes nothing 
and we have the general result 
oH 


Tne OF 


Remark 2: We have used the canonical equations only 
to derive Liouville’s theorem; since then we have not used 
the equations. However, Liouville’s theorem is still true 
if we carry out a transformation whose determinant is 1. 
Thus, if we let 


N 
; = Bix(Q) Px » where [Bul=1, 
=1 


then all of our results are still correct if (7, q) is substituted 
for (p,q). In particular, we have 


oH 
te Or; 


=06,. 


c. Entropy 


We consider a system (#, V) which consists of two ener- 
getically independent subsystems (#,, Vi) and (£;, V,). 
Then, the Hamiltonian separates in the same way: 


H = H, (a) + H,(2®). 


Let W(E,, V:; E., V2) be the probability that system lis on 
the energy shell of the phase space corresponding to(E,, V,), 
and that system 2 is on the energy shell of the phase space 
corresponding to (#,, V,). Then, the probability W is pro- 
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portional to the product of the shell volumes 


w,(H,) dE, w(H,) dz 
W (By, V3 By, Va) aR, = 
_— @1(E,) @2(H,) 
7 w(B) a 
where 
E=H#H,+4+ #£, 


and the denominator assures correct normalization. For 
fixed E, the most probable value of H#, is determined 
from dW=0. Also, for fixed H, d#,——dE,. Therefore, 


puee: 3 dite 5 or ie ie 
* de, * dz, : w,4aE, ow, dB, 
If we define 
01,.= loga,., 
and 
1 doy, 
a 
then 
6, = 6, 


That is, 6 has the properties of a temperature, o has the 
properties of an entropy, and we are justified in writing 


ko = S+ constant 
and 
6=kT., 


By going to a special case, for example, that of an ideal 
gas, we can determine that k is the Boltzmann constant: 


R 
=. 
L 

Remark: Earlier, instead of identifying k as the Boltz- 
mann constant, we took k= 8 and O=kT. However, it 
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can be shown that this is of no consequence, because 


2+—oa = O(logN) 


e—9=0(5). 


10. CANONICAL ENSEMBLE 


and 


As before, we consider a system 1+ 2, consisting of two 
noninteracting parts: 


H= HGH, (p,q)+ 4. (P, 9). 


We let the number of degrees of freedom of the ‘small 
system,’ 1, be much smaller than that of the “arge sys- 
tem,’? 1+2. We now ask about the probability that the 
small system is to be found in the volume element d°”Q, 
of its phase space, independent of the location of the 
second system. Letting 


d?*:2 = w,(H,) dH, , 


£,<H,<E,+02, 
we have 
ad?" 2, 0,(E,) dk, 


a a fd? Q, w,(E2) dz,’ 
where the integration is to include the values of E, = E—E, 
for which F,< FE: 
nie d?*: 0, 02(H — Ey) = a’: 2, w,(H — E,)/w,(E) 
~ [d?:Q,0,(B— By) [d?2,0,(E — E,)/w,(B) * 
LXE E,<z 
Defining w,(F)=exp[o.(#)], we can write 
d?”: Q, exp[o,(# — E,) — 9,(E)] 


oe [a Q, exp[o,(Z — #,) — 0,(B)] 


E,<E 


If we now use the assumption V,<WN,, then 


o,(H— H,)— o,(#) = — «#, (1+ O(N, /N,2)] , 
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where 
do,(£) 


a pie eel, 
— a as 


d# 


In this way, we obtain the distribution for the small 
system. Writing H instead of # and dropping the index 1, 
we have 

Aan = exp[—@H(p,q)]__ 
fexp[— «H(p’, q')]4?*.2’ 
This is the probability distribution for a canonical ensemble. 
Thus, a system in contact with a large heat reservoir cor- 
responds to a canonical ensemble. 

If we specialize to the consideration of a molecule in an 
ideal gas, we obtain the Maxwell-Boltzmann distribution 
anew. 

If a canonical system can be separated into two parts, 
such that H=H,+H,, then 


C= C103 s 


which means that the subsystems are distributed canon- 
ically also. This, together with e=o(H), is characteristic 
of the canonical ensemble. 


11. THERMODYNAMIC FUNCTIONS OF STATE 
Following Gibbs, the free energy F' is defined by 


exp[— aF'] = jexp[— «H]d”Q=Z, 


or 
FP=—kTlogZ. 


The function Z is called the partition function. We can see 
that this is a suitable definition of the free energy by 
showing that the following two equations, which are re- 
quired by thermodynamics, are satisfied: 


) oF 
s=-G, earr(), 
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Differentiating 
Jexptear—myarra =e] 
with respect to «, and using 0H/O#=0, we obtain 


oF 
[(«s+?-#) exp[a(F— H)|d”Q=0. 
From this, remembering the relation 9@=exp[«(F—H)], 
we obtain 

oF 


a+F-—-H=0 
0a 
or 
OF 
—PfaatPr-E=0, 


where H=E. We see that the definition of entropy which 


must be used is 
OF 


saa Ty 


Let us compare this entropy, whose definition is based 
on the canonical ensemble, with the entropy defined by 
means of the microcanonical ensemble: 


exp [— af] =|exp[— aH]d?*72Q. 


First, we integrate over the energy shell. For this purpose, 
we write 


exp[— af] = | exp [— «H]w(H)dH =fexp [o(H) — eH]dH. 


The integral can be evaluated by the method of steepest 
descent. The integrand has its maximum at H=H: 


(az) =< (EH = the most probable energy) . 
dH) z 
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Therefore, 

oH) ee = ony aE +5 ( 
which leads to 


do 
dH? 


) a Ey? + 0[(H—B)'], 


E d2a 


exp[— af ]~ exp [o(£) — «E]] exp 5 dB (Hq — By dH , 


where the approximation is justified by the fact that the 
exponential falls off extremely sharply. (Of course, we must 
require d?¢/dH?<0, which is in agreement with thermo- 
dynamics.) Thus, 


exp|— «F'] = exp[o(E) — «#] V aa 


or 
27 
—aF=o—aHh ] ———_—_ }], 
oP=o—aB+ poe (ecame) 
As a result, 
k 27 
can — Sean —log | ————. 
San = Sninon + 5108 ( aeam) 
or 


= Scaa — onan — (log V) . 


The canonical entropy, as we have defined it, differs from 
the microcanonical entropy only in an unessential way. 


a. Energy distribution 
By definition, 
fw—m expilia(f'— H)| d*7 G6" 
Differentiating with respect to «, we obtain 


o# oF 
set |E-H) («5+ -#) exp[a(F—H)]a*?7Q=0, 


or 


5 + | BH) expla Hy] a" = 0. 
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Thus, to order O(H?/N), the mean-square fluctuation of 
the energy about its most probable value is 


(H—H)? = kT? — = kT C,, 


where C,, is the specific heat. 


b. Equipartition and virial theorems 


As with the microcanonical ensemble, it can be shown that 


on 1 


Oe ae ik y 
where =p or q. This equation is an expression of the 
equipartition and virial theorems (the latter, of course, 
only when the interaction with the walls is included 
in H [A-3]). For example, we obtain by partial integration 


0H | p,(0H/Op,) exp[— «H] 4272 
Pian, fexp[— aH ]a**2 


= (1/a)[p. (0/Op,) (exp [— aH]) avo 
= fexp[— oH | d#¥Q 


RIB 


c. Processes with a variable external parameter 


Let H=H(p, q; a) depend on a force parameter a. If we 
‘ take the logarithm of 


exp [— aF] =| exp[— «H]d?7Q 


and differentiate with respect to a, we obtain 


oF [(0H/0a)>2exp[— aH] 472 
éa  ~—sfexp[— aH] d*2 


or 
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The work done by the system is 


Asa =o (5a) oe 
0a} pa 
hence, 
oH 
es 


For example: a= V, A=>p, etc. 


ad. Entropy 
If we write S=ka(H—F) and E =H, then 


28) _ pq (°2_ (H 
Oa)a ‘ oa oa mel 
In addition, we can vary «: 
~8 = aH + tog {exp[— aH) dQ , 


— da|H exp[— «H]d?*Q 
f exp(— «H]d2*7Q 


~as = deH+ 


+adH — «| (0H /a) exp([— «H]d24Q 


| Chap. 2 


f exp[— aH]d2*Q 


Since the first two terms cancel, we obtain 


a oH 
7% 08 =) jaz — (=). aa| 


or 


1 = (0H 
7,05 = a jai — Gly aa] : 


a. 


The first term on the right is the total change in the in- 
ternal energy, and the second is the work done by the 
system (the change in the energy which results from a 


change in a). 


Because of this result, the condition for adiabatic changes, 
dS=0, shows that the internal energy changes only when 


work is done. 
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12. GENERAL DENSITY DISTRIBUTIONS 


Let @ be a normalized density in phase space: 


feas@=n. 


Then, following Gibbs, a generalized ¥ function (see Sec. 6) 
ean be defined as 


# =e logo d?72 = logo. 


For the canonical distribution, 


oe = exp[a(F—H)], 
and 
8 


#=o(F—H)=a(F—-H)=—-—, 


which justifies our definition. Thus, we can ascribe an 
entropy to every statistical ensemble, even to non- 
Stationary ones. 


a. Theorem: # is time-independent 


From 


eo ee ey (ieo 
Ot =-> (52 Ope Ope Ox) ’ ot : 


it follows that 


OF _fiogo 28 arn0-—— ¥ {tone (22 22 — 22 28) ane 
nl ead ad 2 | loge Og: Op, Op 0 

) oH ) al 
= a = Mae ] or. cane a27Q 


g é 2N 
=-3 {x -[(etoge- os|- sa |(eloee — 0) sha Q. 
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This expression vanishes, as can be shown by partial inte- 
gration, if the forces exerted by the wall are taken into 
account in H. 


b. Lemma 
Let 
L(«, y) = «(loge —logy)—2+y. 


If « and y are positive, then L(#,y)=0 when z=y, and 
L(2,y)>0 when 7¥y. 


Proof: 
oe = log x — lo 
ae s BY, 
ob 1 
ae > 
ae 7a when 720. 


c. Theorem concerning the average 


Let Z be a region in phase space, and let a ‘‘coarse den- 
sity”? be defined by 


pal 
= Gleara, 


Zz 


where Of ial, 
Z 
= [eloggaa, H = Qo logs. 
z 


Since 


ec a_i a) dQ = 0 ’ 
z 
therefore, 


#-# =[o log 0 d242 — 26 log 6 =|Z16, 6)d7Q2>0. 
Zz Zz 


Thus, 
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except when 9 =4, in which case #=#H. It can be hoped 
that the #, which is formed from the coarse density, de- 
creases. This can be proved only in quantum mechanics. 


d. Separation theorem 


Let the total system consist of two parts, 1 and 2, such 
that d”/2=d"0,d?"0Q,. Let the density contain cor- 
relations given by 0(%, #,). By integration we construct 


01 =o d*., and Qe =e ert, 


Since 
Jeroen =fe a0, | ode0, =a" 
therefore, 
[ee — 0:02,)d72=0. 
Let 
HA, = 1 0, log 0, MO Fe 
HK, =o log 9, d2¥22, , 
# =[ologe a 
Assertion: 


KH —H,—H»>0, 


where the equality holds only when @ is a product (no cor- 
relations). 


Proof: 
KH —H,—H, =| efloge — log 0,0.) d?*22 


= | L(@, e102.) d*2>0. 
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e. Combination theorem 


We consider a collection of ensembles numbered by the 
index 7, where 


Joao — 1 ° 


Let o= Dc¢.o;, where }e;=1 and ¢;>0. Then 
a c.f log 0,02" —fe logo d2*2 


= ¥e,| es(loge,— loge) a2 


— =e L(o:, 0) ad? Q>0. 


Thus, 
ASE SU, 


where the equality holds only if 0,—o, for all 7 and k. 
When the ensembles are combined, # decreases. 


f. Minimum property 


For fixed H= | Hod?" Q, # is a minimum for the canonical 
ensemble. 


~ 


Proof: Let 0,=Aexp[—e«H]. Then 
fe log e,d?*2 = log A — [He d?7Q2 = log A— ak. 
Fixing H implies the subsidiary condition 
{ o log 9,d*#*2 =e log 0, 42°22. 
Then, 


AH =fe logo d24,Q2 — fe log 0, d272 =| Lie, Oa OS. 
Q.E.D. 


It should be noted that, instead of the energy, the average 
value of any quantity could have been specified. 
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13. APPLICATIONS 


a. Diatomice molecules (dumbbells) 


Bua 5 (68+ psin®s), 


CM ein : 
P = ; = Af 
6 a6 ’ 

OF kin ae 
Po= caer A sin?6. 


If the potential energy is zero, 


eis a a! 
ae (vt-+ 9% a) 
The volume element of the phase space corresponding to 
6 and @ is 
dQ = dp,dp, d6dp. 


If we define 2,=p,, m=p,/sin@ (angular velocities about 
9, and ~,; see Fig. 13.1), and w=cos6, then 


d‘2 = dz,dz,dudy 


~, parallel z axis 

9, in (x, y)-plane 
perpendicular to 
molecule 
axis MM 

gy = angle perpen- 
dicular to @, 

# = angle perpen- 
dicular to dy 


Figure 13.1 
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and 
1 
H= oA {i + m3} + EF yor(u, y) i: 
_ Because 
oH oH 
— =kT 
Pe Op Po op ’ 
we have 
oe pe ee 
we 
or 
W Opes = ee 


This model is obeyed very well for diatomic molecules. If, 
further, we take into account the translational energy, 
Ev = 3kT, then 


By, = §kT 
It follows for the specific heats that 
Cot 
C=, Cp sek, Wand yeq ayo. 


At low temperatures the rotational degrees of freedom of 
the hydrogen molecule are “frozen in.’’ This is connected 
with the fact that classical statistics arbitrarily counts 
certain degrees of freedom and does not count others. For 
example, here we have not counted the degree of freedom 
corresponding to rotation about the symmetry axis. 


b. Rotating rigid bodies (general case) 


We introduce the components of angular velocity along 
the principal axes, expressed by Euler angles: 
u, = sind siny+ 6 cosy, 
u, = y sin cosy — 6 siny , 
Us = ycosd+ yp. 
In terms of the principal moments of inertia, A,, <A,, 
and A,, we have 


Fyn = $(Ayui + A,ug + A; 43) : 
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Hence 

Po = A, cospu, — A, sinpu, , 

Po = A, Sin sinyu, + A, sin# cos yu, + A; cos6 us , 

Dy =-AzUs 9 


and 


O(Pis Pas Ps) _ 4 4 Asin d. 
O(Uy y Us y Us) 


Then, 
dW = constant x exp - wr ds2 


= constant x exp |— (A, ui + A,uj + A,uz + 2 Boot) 


st 
X du, du,du; sin dédydy, 
and 
4A, ui = 4A, u} = 4 A,u2 = 4 kT. 


For polyatomic gases, we have 
rotational energy = 3kT, 


translational energy = kT, 


and 
E=8kT. 
Thus, 
C, 
Goi, “C,—A4n, and le 


e. Vibrating dumbbell 


Introducing the reduced mass 


#=Mm|(M+m), 
we have 


2 
Ge ep wl te seb Bese « 
Qu 2 


For each vibrational degree of freedom, H=k&T. This con- 
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tribution to H is in addition to the contributions from 
other degrees of freedom. 

It is noteworthy that this energy does not in any way 
depend on the coupling strength (does not depend on @,). 
Therefore, for w,—co (rigid dumbbell), the contribution 
certainly does not go to zero. However, it is known experi- 
mentally that this degree of freedom generally does not 
enter. 


d. Solids 


We imagine that a solid consists of atoms which can 
oscillate about their equilibrium positions. For small oscil- 
lations the force will depend linearly on the displacement 
and, therefore, the potential energy will be a quadratic func- 
tion of the displacement. Therefore, according to Euler’s 
theorem, 


OF sot 
3 Ik Og 7 cer 
and 
Eyot = kT -3N. 
Also, 


LL =r Sy, 
which leads to 
- E=3RT. 


Thus, C,=3R—~6 cal/°C. This is the law of Dulong and 
Petit. 

Deviations should be expected at high temperatures, since 
the displacements at high temperatures will not be small. 
However, the law is well satisfied at high temperatures, 
whereas, since C,—>0 as T—0, it is not satisfied at all at 
low temperatures. The explanation of this requires quan- ' 
tum statistics (see p. 82). 


e. Langevin’s theory of paramagnetism 


We consider an ideal gas in a magnetic field H, each of 
whose molecules has a permanent magnetic moment p. 
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(Exactly the same considerations are valid for dilute so- 
lutions.) We need not discuss the kinetic energy, which 
only depends on the z’s; 


Egor = — wH cos 6 , 


where 6 is the angle between yw and the magnetic field H. 
Therefore, 


W(6)d6 = constant x exp |- rd sin6 dé, 


which leads to 


feos exp [— (Eyo./kT)] sin 6 a6 
M = 0086 = wp; ; 

fexp[— (Pror/kT)] sin 6 46 

a 
From the partition function, 

F r po . 
exp |- al = [exp E cos | sinf dé, 
0 


it follows that 


kT wb _ #H)\_ , sinha 
= Fa (e [ee| | male? 2 


Therefore 
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Limiting cases: 


(a) a>1: M=p (saturation) . 
a 
(b) a@<1: F~—iES +... 
ae 
tt eee 
The susceptibility per mole is 
_lw,_ 1 (uly 
A~ShE 3 HP 


This is Curie’s law. 


| Chap. 2 


Remark: It is essential to ascribe a permanent magnetic 
moment uw to the atom without inquiring about its origin; 
only quantum theory can answer the question of the origin 
of the permanent magnetic moment. As a matter of fact, 
if the magnetic moments are considered to arise from ele- 
mentary amperian currents, and if the elementary parti- 
cles (nuclei and electrons) are treated statistically, then 
the absurd result follows that paramagnetism does not exist. 

For charged particles in a magnetic field H, we have 


where 
H= curla 
and 
0A 
ia OF 


The absence of paramagnetism can be verified as follows: 
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Hence, 
pO ee OL sO € f 
ma, = ; Dull (= = a) > be ) 
me € : 
mg = — Hxq, 
mv =—vxH (Lorentz force). 
If we define 
m= p—-A — mq, 
then 
d*pd*q = d®xz dq , 
and 


alt 
d*v W = constant x exp |- a — nan 
2m 
m 
= constant x exp | as | dv, 


where W is the probability density for finding velocity v. 
This holds whether or not we assume a potential energy 
between the particles. Thus, we always obtain the Max- 
weil distribution and, in particular, there is no magnetic 
moment. The reason for not having a magnetic moment 
is illustrated schematically in Fig. 13.2. Those particles, 


Figure 13.2 


indicated by a dashed line, which are reflected by the 
walls, give rise to an oppositely directed current that can- 
cels the moment produced by the other particles. Para- 
magnetism is a pure quantum effect. 
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f. The van der Waals gas 


For the sake of simplicity, we restrict our attention to 
a monatomic gas. Let the potential that describes the 
internal force between the ith and kth molecule be U(r), 
where U(co)=0. Then (NY = number of particles) 


Exot = = > Eh 
i<k 
1 


N 


es aa Enor 


M Kor 


exp [—aF] a lfeel aon P| 


<exp[— aB pot] d8p, ... d®py dq, ... d°qy . 


If we ask just for the velocity distribution, then we again 
obtain the Maxwell distribution 


exp[— a«F] = (2amie [exp [— ogee Gges od ae. 


We cannot evaluate this integral exactly. However, we can 
make an expansion that is essentially an expansion in de- 
creasing powers of V. Let 


W,x(1 ix) = exp {— aU (r;x)] == 1 3 W,.(00) => 0 5 
We expand 
exp [— aye] = [[ (1+ Wi) 


i<k 


=14+ )Wa+ ~ Wa Wiw+.. . 


i<k (§<k) 4 (i' <k’) 


The linear terms contribute only when two atoms are near 
one another; the higher-order terms contribute only when 
more than two atoms, or two or more pairs, “collide.” 
For example: 


W,.W,,;: (123) together; 
W,,W.,: (12) and (34) together at the same time. 
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Integrating term by term, and introducing center-of-mass 
coordinates, we obtain 


iz (12) d°q, ... d°qy = a dq, dq, 


co 


0 


| Wi; d3y, = Veale . 
Therefore, 


fexp [— CB 0+ | d*q, vee d*qy 


Analogously, 


© 


Vay (; | ocr? W(r) dr + O(V¥-2) 
0 


~ V4 a Azer? W(r) dr + O(V*-2) ; 
. | 


exp [— af] = (2amkT)34/2 VX 


x1 pe fomm r)dr+0O Fal : 


i = 4 N log (2amkT) + N log V 


kT 


ee 
+ ry fore )4ar2dr + O (F:) : 
(1/V)> A(P)/kT 
This gives per mole, with p= ‘wae 
_ kT am ) 
a : 
This is the same result obtained with the virial theorem in 
the kinetic theory of gases [A-3].? 


1 References on higher approximations: H. D. URSELL, Proc. Cambridge. 
Phil. Soc. 28, 685 (1927); G. E. UHLENBECK and E, BETH, Physica 3, 729 (1963) 
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Figure 13.3 


If U(r) is as sketched in Fig. 13.3, then, as in kinetic 
theory [A-3], one finds 


A = Fie, pe 


2 
a= — = [ ert a: 


RT(, b\ a : 
D a (1 7) V2 (van der Waals). 
14. GRAND CANONICAL ENSEMBLE 

a. Homogeneous systems 


In a homogeneous system of volume V and WN identical 
molecules, we separate a subsystem of volume V, contain- 
ing NV, molecules by means of an imaginary wall (Fig. 14.1): 


V=V,+ V,, NV == Nee, 


Figure 14.1 
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We want to determine the probability of finding N, mol- 
ecules with coordinates q’,...,q" in V,. First we ask for 
the probability of finding WN, specific molecules in V,. We 
begin with 
W(M,, Vi, x) = exp[+ 2F(V, N, «)] 

x [ fexp air, on, a pm, nn, PIA, LHD, 

V, Ve 
where 

exp[— aF(V, N, «)]=|exp[— «H(p, q)]d2472. 

Vv 


If we assume that V, is large, in order that surface effects 
can be neglected, then we can neglect the interaction be- 
tween the molecules in 1 and those in 2: 


Pap, .... eer, ...0%), 
| ve lie wHy(p', ..., g%1)]d2™Q, 


iV; 


xfexp [— «H,(p™*1, ..., g%)]d?%Q, . 


Ve 
Therefore, for N, specific molecules in V,, 
W(N,, Vi, «) =exp|a{7(V, N) —F,(Vi, N,) — F3(V2, N,)}] ) 


where F, F,, and F, are the same function [A-4]. 
Now we return to the original task of determining the 
probability of finding any N, molecules in V,: 


W*(N,, Vi; &) 
N! 
alg: iM! ae ee (VN) F(,N,) FAV, N2)}) - 
We introduce a new function #* by 


a exp[—aF]=exp[—oF*], F*=F+kT logN!. 
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Then, 

W*(N,, Vi)= exp[a{F*(V, Ny) Sa N,) —Fi(V., N —N,}). 
The condition for the most probable N, is 


log W* = maximum 


or 
oft OF; 


ane ONG 


If we define u=—0F*/ON, then the condition is 


Ma(Vi, N;) = H2( Ve, Ne), 
where 
N=N,+N,. 


Fluctuations. We begin with the Taylor expansion of 
log W* about the most probable value: 


a OFT « = 


log W* = log W,— 2 ON? (AN,)? — 


By using AN,=—AN,, the expansion can be written as 


OFT | OF! 
ONi ON: 


log We= log Ws — 2 / amped. 

The higher-order terms are small if the average values 
of NV, and N, are large [A-5] (which was already assumed 
in neglecting the surface effects). Therefore, as before 
(compare p. 23 with the ideal gas result below, for VY, < NV), 


kT 


(m= O2Ft/ON? + OOF /ON2 


Application to an ideal gas: 
—akF = NlogV+ Nf(a), 
al* —logN!— N log V— N f(a) 
= NlogN—N—Niog V—Njf(a) ; 
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hence, 
au = log N — log V— f(a) . 


Since 4,—=,, therefore, V,/V,=N.,/V,. From 
ot Nt = Tan 


- it then follows that 


TS a a ee 
ee ennagee 


In this case, F* is a homogeneous function of the first 
degree in V and N. It can be shown that this is also true 
in general: 


PCN ANG gir my IN, o) 
Using Euler’s relation for homogeneous functions, we obtain 


ore _OF* 
es ea, ee 
ae oe 


hence, with oF*/oV = — p, 


Nu=F*+pvV 
and 
F*+pV 


which is known from thermodynamics. Also, uw is a homo- 
geneous function of degree zero in V and XN: 


BAY, AN, «) = w(V, N, x). 


Euler’s relation implies 


OH yy SH 
aN ; 
or 
o2r* O2F* O2F* V Op 


me 
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Since p also is homogeneous of degree zero, therefore, 


Op Op _ 
N aN 4- v ay = 0, 
which leads to 

on V\2 op 

on2 ~~ s\N) OV" 

In order that W, really be a maximum, we must first have 
O2F*/ON*>0, 


which means 


This is a stability condition. Second, because of homoge- 
neity, we have 


(ee a a ee 


‘OV, °OoV, OV’ 
or 
kT 
AW,)?= — 
(OMY = — (VIN OplOV) VV, + 1/V;) 
With 
Bee V op 1 op 
n= 7 and WON Nn? 
this becomes 
a 


Limiting case: 


Vie, = Ve Nv< N, 


—_—— eas tl 
AN,)? = ———— N,. 
(OM! = rem 


b. Generalization to mixtures 


Let the index label the components of the mixture. Let 
there be WV i molecules of component 7 in the total vol- 
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ume JV, and let there be V, molecules of component 7 in 
the partial volume V which is separated by an imaginary 
wall. The associated probability is 


py *(N,; Ng, --., V0) 
it = 
2 Nye ecu, ON, Nast -» V) 
v 
— F*(N,— M,,..., V—V) —H(p', ..., @)]} a2, 
where 
F* = F+kT > logy,! . 
k 
For the most probable distribution, 
Ae = — oF* 
ee 
In the following, we always assume V< V. Insofar as 
we restrict attention to the linear terms in N,, we have 
_—- 
NENG 
x fexp [«{2+ uN, +.N,+...—H(p',..., g)}] dQ, 


Vv 


W*(N,, N,, eoeg V;, ) aad 


where we have put 
oF* 


ay ae 


Q=+ 


For equilibrium, there is the relation 


2K 
2 —p=—p or Q=—pVv. 
Consider 2 as a function of the px, Q(p1, Me, «3 V, T). 
For fixed 7, because of the homogeneity, there is the 
relation 
oFr* 


r= atl Sp: 
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Therefore, 


Q=F*—~ Nim. 
12 


From this formula we can derive 


(=2),= (er), 


The quantity 2 is a new thermodynamic potential. 


Averages and fluctuations: 


DECC Peon fe aa 
Ne 


can also be considered to be the definition of Q: 


Chap. 2 


exp[— 20] = 5 DE ta” a H}}. 


Differentiating with respect to w,, we obtain 


ex) 


ae lhe 


Ny 
Therefore, as before, 
a? == 
Ou: 


Differentiating again, we obtain 
=5 1 Ea 


From this follows 


Of: Ope 


02.2 


AN‘ANt = —kT 
OM OMy | 


veal ge 
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This result can also be obtained with the help of the fol- 
lowing form of the probability, which follows for V<V, 


w* = we exp|—5 AN. | 


5 Sayan Saran _ 
and with the definition 
ANNA =|. ; faax) ANA 


The calculation of the integral is easy if the following 
lemma is used: 


Lemma: Tf 
a — i> 9x2‘ X* , 
i,k 
then 
a ie .{daxtx* exp aie) 1... 
[-. . [dx exp [— xf(x)] 7 
where 


297 Gur = Ob. 
This can be proved as follows. First we show that 
ae 
one a 
This follows from the relation 


ae - —* 


Cn 


by partial integration. Therefore, 


' 
er ; re 
ee 
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which leads to 
=> eee 
gigt=—g*, Q.E.D. 
4 


The above result follows with the identification 


_ OF  m& 82 aM 
de O8N,ON, ON,’  duidm on, 


Gre I 


By differentiating with respect to «, combined fluctuations 
of particles and energy can be calculated. 


Chapter 3. Brownian Motion 


15. INTRODUCTION 


The only observable quantity is the displacement of a 
Brownian particle, 


A(t) = x(t) — 2(0); 


its velocity is not observable. As is always the case with 
irregular motions, the mean square of the displacement 
will be a linear function of time: 


A(t)? = 2Dt. 


It turns out that the distribution of 4 itself is Gaussian: 


1 Ae? 
W(A)dA = ———e |- il a4. 
eo en | ah 


In the same way, rotations about an axis can be considered: 


A) = at. 


It should be noted that Brownian motion cannot be used 
as perpetual motion of the second kind; this is because 
the phenomenon is exhibited by all physical devices. 

In this chapter three theories of Brownian motion will 
be treated. They are the theories of 
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1. Langevin (with the virial theorem), 

2. H. A. Lorentz, and 

3. Einstein (comparison with diffusion). 

It suffices to consider the linear case, because the motions 
in different directions are independent of one another. 


16. LANGEVIN’S THEORY 


We separate the force on a particle into an ordered part 
and a disordered part. The ordered part is a resistive 
force equal to —Wz. The disordered part arises from col- 
lisions between the molecules, and it will be designated by X: 


Because of the disorder, X=0 for fixed x or &. Now, 


mor = — Wae-+aX, 


m < (a) — me = — oe ai (wv?) + aX, 
If we take the statistical average of this equation over 
many particles, then one of the terms, xX, will be zero; 
this is because the collisions at a fixed position occur in 
a completely random manner. Furthermore, the operations 
represented by d/dt and the overbar (statistical average) 
commute: 


RCE wo d PSE 
— a —, 2 = — — 
ma (xxv) — max Ww a ($@?). 
According to statistical mechanics, 


me —2E. = kT. 
Thus, 


ry = (ip) +495 (a2) = kT. 


Let us choose the origin such that 7(0)—0. Then the so- 


Sec. 16 [| LANGEVIN’S THEORY 65 
lution is 
= (x2) = ae + Cexp |- 4] : 
Since W/m =~10~ sec-', the stationary state is reached very 
quickly, and we obtain 
ey. 
W 
17. THEORY OF H. A. LORENTZ 
Let v=; then 
mv=—WrtX. 
Integrating this equation from 0 to t, we obtain 


mM(V:— V) =— Wut+G,, 


Ge =|xar : 
0 


For sufficiently small ¢, using the fact that —Wv is random, 
we obtain 


where 


Vv, = VU (1-Te)+o@, 
m 
( —— (120 a) 42 So, (1- Te) 4S Se esse.) 


If we take the statistical average of this equation and 
note that v2 must equal v? and that G,v,=0 (as before), 
then we obtain 


Gi = 2Wmvst = 2WkTt. 
This is the mean square of the momentum transferred to 


a particle in the (short) time t. 
If we define B by f=1—Wt/m, then 


1 
v= Brot — G.. 
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Consider a time interval of length nt. and let t, and @, 
be values at time ki, k=—1,.2,....”. Then 


iL 
Vv, = Br eg 


3 1 
v, = Br, Tas = Soe (8G, + G,), 


b! 
= By) + = (Br-G, + B*-*G@, +... +G,) . 


Letting v= G,/m, this can be written as 


1 
Cn =a + Le a ( 
m 


The quantity Av can be written in the following way: 
Ax = tgp +0, +... + Va-1) 
wai a (1+ B+...4+ 6") 
+ = G14 B+ .04 8) $e t oOo 


Go1—p? t 1-9 t —_~p 
— Bp ( —_ — St -- s. 
mip am 1—p m 1—~p 


By using the relations 
Ga%M=0 and GiGy=danG® (nO), 


the mean square of Ax can be written as: 


mS Cp) rw) 


yal 
gil ak o.oo, 1—”" 3 
G2 ————__ (1— 26"+ 3s") 2- —_ @& -) 
@ 2a “Wr ilies m: ad Geer 
dll B ? Qa _ L 


(1—B)? "m= (1— 8p)? 


1p _ se 
x(n—1 28 ; 3 _p! — }. 
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We now assume n>>1. Then the term proportional to n 
is the most important. With 1—B=(W/m)t, that term 
becomes 


— 
(W/m) OW 


2 
(Any = — ower 
m 


This result is the same as that predicted by Langevin’s 
theory. 


18. EINSTEIN’S THEORY 


This theory, which was put forth in 1905, is the oldest. 
According to statistical mechanics, the suspended particles 
exert a pressure 


p— wel , 


where » is the number of particles per cm’. The diffusion 
coefficient D is defined phenomenologically by 


i=——Dgradn, 


where i is the current density. In one dimension, this 
equation is 


If an external force K and a resistance —Wv act on a 
particle, then, in the stationary case, 


C= kW, 
which leads to 
[an ais 
Ww 
In a stationary state in which there is current flow aS a 
result of a pressure gradient, the force per unit volume 
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must maintain equilibrium: 


dp dn 
te 
Therefore, 
kT dn kT 
Sede ee 


Now, (Az)? will be related to D. There is a probability 
density g(s) that a particle suffers a displacement s in 
time rt as a result of Brownian motion. Let the number 
of particles between s and s+ds be 


dn = ny(s) ds. 
Then 


+00 


nav, t-+7) dz = dn n(x — s, t)g(s)ds. 


For tv small, 


+c +0 
ra) 
n(x, t) + = ae fois ds — = fst) ds 


If we assume that the higher-order terms are small to a 
higher order in t, then 


Because of 


the diffusion equation becomes 


an_ etn 
ota? 
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Comparing the two equations for On/dt, we obtain 


Brownian rotation 


This can be treated in a completely analogous way: 


: g? = — a? = 4kT 
This leads to 
—= 2kT 
Pt = Ww io: 
With 
t 
fi — Oar, 


0 
the result is, as in the theory of Lorentz, 


I? =2kTWt. 


Example: Current fluctuations in a closed circuit without 
electromotive force. 


t 
With G,= | Xdt, the calculation proceeds as in the theory 
0 


of Lorentz: 
Gi=2RKTt, 4Li? =4kT, 


Chapter 4. Quantum Statistics 


19. THEORY OF BLACKBODY RADIATION 


On purely thermodynamic grounds, it is possible to derive 
the following two laws:? 
1. The Stefan-Boltzmann law for the total energy density: 
E 


u=Zaal'; 


2. Wien’s displacement law for the spectral energy density: 


=e (=) : 


For calculating F(v/T) we have to use statistical methods, 
of which there are two: the method of oscillators and the 
method of normal modes. 


a; Method of oscillators 


By an oscillator we shall mean a mechanical system 
(a harmonic oscillator) which can emit and absorb a spec- 
ified frequency. On thermodynamic grounds, the equilib- 
rium radiation must not depend on the details of the 
structure of these oscillators. Therefore, we may immedi- 


1 See W. PAULI, Lectures in Physics: Thermodynamics and the Kinetic Theory 
of Gases (M.I.T. Press, Cambridge, Mass., 1972). 
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ately assume that the oscillators are harmonic oscillators: 
= ge + Le @? 2 
2m 2 q 


It can be shown that such an oscillator in contact with 
radiation has an average energy 


= c 
"= Baye 

Furthermore, according to classical statistical mechanics, 
the relation 


ee 


should be valid. Substituting this into the previous for- 
mula gives 


Srv? 
Ov = ce kT ’ 


which is the Rayleigh-Jeans formula. It can be seen imme- 


diately that this formula is incorrect because fi 0, dy = oo 
(the ultraviolet catastrophe). 

In order to overcome this dilemma, Planck made the 
following hypotheses: 


1. The energy of an oscillator can only assume discrete 
values: 


HE, = Ey+ne(v). 


Also, the energy of radiation is absorbed and emitted 
only in small quantities n’e. 
2. The relation 


is kept. As a consequence, the statistical mechanical 
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partition function is now a sum instead of an integral: 
Z => 9,exp[— LE, /kT] . 
n 


Here the third hypothesis is natural. 
3. The degeneracy of state n is g,. (That is, state is 
gn-fold degenerate.) 


Corresponding to these hypotheses, the microcanonical 

ensemble is to be defined by 

a E<#H,<E+dE, 

constant x W, = ? 

0, otherwise . 
Further, if the total number of states in the energy shell 
(taking the weight of each state into account) is U, then 
the entropy is 

S= koe VU. 

With these new concepts, all further formulas of classical 
statistical mechanics are still valid. For the case of the 
Planck oscillator, remembering that its states cannot be 
degenerate,? we obtain 

Z = exp[— E,/kT] > exp[— ne/kT] 


1 
— exp[— e/kT] 


= exp(— #,/kT| i 
per oscillator, or 
P= pig r= kT log(1 ~ exp[— e/kT}) . 
Again taking 


Ck, 
we have 
= 6 6) 
B= (al) =— x logZ 
or 
a E 
°" exp[ae]— 1° 


* See W. Paul, Lectures in Physics: Wave Mechanics (M.I.T. Press, Cam- 
bridge, Mass., 1972). 


Sec. 19 | THEORY OF BLACKBODY RADIATION 73 
Using hypothesis 2, we obtain 


- _ Si E 
eS Cepia |— 1° 


This has the form of the Wien law, 


V 
o,= VF (7) 4 


the validity of which is assured by its thermodynamic 
derivation only if 
e=—hy, 
where h is a new universal constant, the quantum of action. 
Finally, 
_ 8ah vy 
ore expl[hy/kT]—1° 


which is Planck’s law of radiation, and 


hy 
exp[hy/kT]—1— 


E=E,+ 
Remark: In contrast with the classical case, the dimen- 


sions of the expressions of quantum statistical mechanics 
are correct. Let us consider the limit hy <kT: 


hy\?2 | 


so that 
y 1 hy 
F = E,—kT log (2 al 
hy 1 
= #,+ kT log i 3 hy, 
and 


In the classical case, we found 


E=kT and een (108% aa constant). 
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We now show that, asymptotically for large T, 

Zauant~ Zeciaee/h : 
Per oscillator, 


Zeuan. = 2 eXP[— nhv/kT] exp[— E,/kT], 
Zeiss =| Ap dq exp[— H/kT]. 


First we integrate over the energy shell H<H< H+dlL: 


dQ 
oab=|[apaq— Fe ae 
E<H<E+{dE 
Using 
mes ue 242 — 
alae aie oe 
we obtain 
ena ae! 
A | 
@ m @ v 
Therefore, 
10 = ap == 


and 
hese =[exp [— E/kT}] ae 


Zavant = > exp[— nhy/kT] exp[— E,/kT] 


=} exp[— vhr/kT | dx = | exp[— E/kT}] — ; 


Letting U be the number of states in the energy shell, 
we have 


U=wdE/h, 


which says that, for the linear harmonic oscillator, the 
density of states in phase space is 1/h. This is a general 
result: For a system with f degrees of freedom, the den- 
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sity of states in phase space, asymptotically for large 7, is 


ome tS 


Therefore, all classical partition functions must be mul- 
tiplied by h7~’. 

For radiation the zero point energy EH, is not impor- 
tant; for solids, on the other hand, it is important (see the 
next section). We can consider two assumptions: (1) E,=0, 
or (2) H,=thy. With the second assumption, the term 
—thy that appears in some of the preceding formulas is 
cancelled. Quantum mechanics requires assumption 2, 
and that assumption gives the correct result for the experi- 
mentally verified vapor pressure differences between isotopes 
(see the next section). 


b. Method of normal modes 


We consider the normal modes of a cubical cavity of 
length 1 on a side, for waves whose direction cosines are a,. 


We have 

e@+toaztoaf=1 and Ete 
cali i ie 

where s,, S,, and s, are integers greater than or equal to 


zero, and 


a 
Se 
A= Vahp ah + at 
For each set (A, «;) there are two normal modes in the 
radiation cavity that correspond to the two directions of 
polarization. Asymptotically, for A<1, the number of nor- 
mal modes in the wave number interval 


(1/a, 1/a + a(1/A)) 


can be determined by noting that, asymptotically, the 
number of lattice points in a spherical shell in the first 
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octant of s-space equals the volume of the shell: 
(Ag eee p—AnReaR 
=p eg le 
= Bagel) = 4(5) 
4nV _ fl 


where p=polarization factor, R?=)> si, R=2l/A, and 
VV. : 


Note: (1) As long as A*<V, this formula is valid for a 
cavity of any form. (Effects proportional to the surface 
area have been neglected.) (2) The formula is valid for 
an arbitrary dispersion law. 

As a special case, consider light in vacuum: 


2 
No, y+ av) = V = ay, 


Let £, be the average energy of the normal mode in 
thermal] equilibrium. Then, 


=, aiaete 
Vody=Ni,9 bdr) B,, = B, 


If we apply quantum theory to the normal modes, then 
we find 
L,=nb, 
which leads to 
a 
exp[hy/kT|—1 


From this we again obtain Planck’s law, 


= ye 
Or exp[hy/kT]—1’ 
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where 
__ 8h 
— = 4 


i 


Some implications of this law are: 


hy hy : 
i. LPO 1: 9, yr exp |- iz (Wien’s law), 

hy ; 82? (Rayleigh-Jeans’ 
4. ap <li @=—Q bP formula), 


3. The total energy density is 


foe) 


el \* | ode 
See? \ ae 3 
it] 


0 


1 DN hast 
- (=) cara) =r (F) 90 °° 
Therefore, 


870° k4 (Stefan-Boltzmann’s 


ee 7 ne 


From this and the Planck formula, h and k can be deter- 
mined separately: 


h = 6.62 x 10-*’ erg -sec . 


c. Fluctuations 


Substitution of > for | does not alter the statistical 
mechanical formulas: 


—_—— On OH —k 
ai ee (5z),-—? Gabe ere 


In the case of radiation, there is a special simplification. 
This formula is also valid for a small subvolume V, since 
the number of particles does not give rise to a new variable. 
Therefore, the grand canonical ensemble is identical to the 
canonical ensemble. Furthermore, the formula is valid for 
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each frequency interval separately; this is because the fluc- 
tuations, like energy and entropy, are independent: 


— ON ,f 08 \ | =k 
ena (sr), = *(sa7n)) =e 


From Planck’s law it follows that 
hy OS, — hy a Os Ov 
th =] 1 == || == 
k OB, kT og ( 7 m) oe (*+1)- 10 a 
Letting Z,=—0,dyV and S,—s,dvV, we obtain 


hy oS, hy Os, Oy Ov 
¥en woe (2 1) a 


= log (e,+ yr?) — loge,. 


Therefore, 


k 
So = 5 [art y9*) log (Ost 709) — gy log oy— yv* log yr*] 


k Oy Or 
.—blnerin( 9) -ane( 


For the fluctuation we obtain 


or 


~ae... (ee®) O(1/kT)]— 
su = EMD] =a BY 
aU/kT) 1 9 
00> ie 00» [log (e,-+ yv3) — log o,] 
= am ales 
hylotyr o, hy 0,(0,+ yr) 
Therefore, 
Sa hy 
(AB, = (cot +087 wm) avy 
or 


————e c 
(AB,)? = E hy + aa | dy V. 
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This formula is due to Einstein. For h=0 (Rayleigh- 
Jeans), only the second term occurs. It represents inter- 
ference fluctuations.? Only the first term occurs if we 
start from Wien’s law, or if we consider light as a 
classical ideal photon gas; in the latter case, we take a 
photon to be a corpuscle of energy hy which moves with 
the velocity of light. 

In that case, let U equal the number of photons between 
yand v+dy; H,=Uhy. Then, according to p. 23, 


(AE,)? = (hv)? AU? = (hv)?U = WE, , 


which is the first term of the previous formula. 
Another way of writing Einstein’s formula is the following: 


2 
The number of normal modes is = = ie 
The energy per normal mode is e¢,: E,=¢,N. 
Thus 


(AE)? = Ne,(hy + e,) . 


d. Radiation in a dispersive medium: 
Let n = n(A) be the refractive index, v= c/n the phase 
velocity, 
dy c Adn 
= 2S (1 
d(1/A) =( a n 7) 


the group velocity, and A= c/ny the wavelength. Then, the 
number of normal modes between yv and »+dr is 


82 1 

7” y om dy n? 

7 ce =1-+-(A/n)(dn/da) 
Sav? dy 

=a el 


2See H. A. LORENTZ, Les Théories Statistiques en Thermodynamique (B. G. 
Teubner, Leipzig, 1916). 
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Therefore, 
ns ha 


0, = ( @»)seeegnins z ae (A|n)(an/dA) ae ne vU . 


20. THEORY OF SOLIDS 


The model of a solid that we use is that of N coupled 
mass points (atoms). Such a lattice has 3N normal modes. 
It can be shown that these vibrations always have the form 
of waves because of the existence of a translation group 
which is required by the lattice structure. For example, 
for a cubic lattice whose lattice constant is d, we have 


cular 
Un nn, — A exp Ee Fe (11% + My + Ns %y— “| : 


Conditions on the «,, analogous to those in Section 19, are 


27 

ai ot = 8. 

A characteristic difference in this case is that there are 
only a finite number of normal modes. However, for 
l>>A>d, the system is practically a continuum, and its 
normal modes are the elastic vibrations. That is, 


v= Ve; 


A 
where v, is the velocity of the transverse waves (p= 2), 
and v, is the velocity of the longitudinal waves (p=1). If 
a je 1 
o UF OU 
then the number of normal modes in the frequency inter- 
val dy is 
peli, 
N= V = yrdy, 
v3 
For small 4 these considerations are invalid. In that 
case, there is a dispersion, which depends on the specific 
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lattice structure, and which is such that the total number 
of normal modes is 3N. Following Debye, it is a suffi- 
cient approximation to keep our distribution up to some », 
and to cut it off at that frequency: 


implies 
3 _ 3Nv° 
aa 4nV 
Eliminating v, we obtain 
ey an vedy. 
Vp 


Using the quantum theory of the normal modes, together 
with statistical mechanics, we obtain 


hy 


a at: exp[hy/kT]—1° 


Therefore, 


*pD 
a pp ON. hy 
pe a or [> ud (2. + exp[hy/kT] — | ; 
0 


a. High temperatures: hy, [kT <1 
Let us define i), for an arbitrary function f(y), by 


3 
flv) = 5 |? flv) dy 
Then, 
= ~ hy 
B= 3N (B+ Soren) 


For hy, /kT <1, we have 


EB =3NkT+3N (z.,— 5) ; 
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82 
The free energy is 


*D 
9N hy 
He Ee oie ae v8 dy| Boy +h log (1- exp |- iz|)| , 
0 


For hy,/kT<1, we have 


ae 


~ hi 
F* — 3NkT log (=) 4+3N (#.— =) 


b. Low temperatures: hy,/kT >1 
Because of the very rapid exponential decay, only the 
smallest values of » contribute. Therefore, we may inte- 


grate from 0 to co: 
: 3 — htdy ON (kT\*. ot 
ca mee = 

BH : | es [hy/kT]—1 ( h ae 


Vp 
0 


where E,=3N B,, : 


h 
With 0,= “2 and Nk=R, the result is 


Therefore, 


c. General case 
Define 
o fpesde 
F(a) = =|a— . 
0 


For the two extreme values of #, we have 


o<1: E(e) ~1, 


2>1: F(a) = (3/a*)(4/15) . 
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In terms of F, we can write H—E, as 
E-—E,= sere (2) ; 


Neither limiting case depends on our arbitrary definition 
of y,. At low temperatures, the large values of » are mean- 
ingless, and, at high temperatures, we simply have equi- 
partition. It is for this reason that the Debye approxi- 
mation is a good one. 


d. Free energy (neglecting the zero-point energy) 


F* = 3NkT log (1— exp[— h»/kT]) 
*D 


= oNiT + log (1 — exp[— hv/kT}) v2 dy. 
D 


Defining 


—=y and ~ =2%= 


and integrating by parts, we obtain 
JIE sNer = {log (1 — exp[—y]) y*dy = 3NKTG(a). 
0 


The function G(x) is defined as follows: 
3 2 
G(w) = — | log (1— exp[— y]) y*dy 


0 
z 


. es ee a 
ee aR 9) | leg oxnl 9] 


0 


dy 


Blas: yf? 
= log(1— exp [ees ae 


0 


= log(1— exp[—«]) — F(a), 
aG'(e) = F(x) . 


way! 
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It follows that 


peasvera(F) =swea|I0g (1—exp[—@,/T}) == (2) 


For 7<@,,, which means #>1, 


ae SL 
Therefore, 
Te T 
2 ~+E,. 
F 5 Nk os + 
e. Entropy 


In the approximation T<@,, 

An4 ye 
w= R(g)- 

Remarks: (1) The meaning of the asterisk inthe symbol F* 
is that F* is the free energy defined by F +kTlogN! (see 
p. 55). That is, it does not contain the term which arises 
from interchanging the molecules. For this reason, S* is 
zero at zero temperature, corresponding to Nernst’s theorem 
(see p. 93). 

(2) At low temperatures the solid is analogous to the 
blackbody, as ean be seen by replacing the Stefan-Boltz- 
mann constant a by 

32*R/(5V O38). 


f. Zero-point energy, E, 


If we divide the volumes of the elementary cells in phase 
space by h’, then we obtain the following results for an 
ideal, monatomic gas: 4 


Or* 
My = aa = KT (logp, — $ log T— i), 
kTN, 
oe Va 


g 
t= log[(2am)tkth-*] . 


“See W. Paul, Thermodynamics. 
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Let ¢,=£,/N,,,3= zero-point energy per oscillator for a 
solid. Then, for T>> @,, 


SS 


[aoa = 3kT ad 


ap e+ ep 


The energy & consists of two parts: (1) the heat of vapor- 
ization, —A,, and (2) the zero-point energy of an oscil- 
lator, @,,. 

hy) = — Ao ae eo» . 


The vapor pressure follows from 4,,4=,: 


me ie a . 
logp = + §log T + 4+ 3 log ~~ <4 +, (By, — $08). 


For different isotopes, the values of » are different but 
the values of A, are the same: 
A= AG ; Yp|Yp= V M/|M' (since Bee = 1m (2zcv)2 7? = Evyot) . 


Different isotopes would have different vapor pressures at 
high temperatures if @,43h7. Because this is experimen- 
tally not the case, one concludes that [A-6] 


Coy = 3 hy 


for a three-dimensional harmonic oscillator. Again, this is 
in agreement with quantum mechanics. Hence, for 0, < T, 


hy A 
= 5 i eee 
logp = + $logT+7+3 log LORE’ 
For 0,>T7, because 0,/M=O,VM", this result is in- 
correct. 


21. ADIABATIC INVARIANTS 


For adiabatic compression of a radiation cavity, it can 
be shown that [A-7] 
E, E, 


Y y 
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Ehrenfest made the general postulate that the quantum 
state of a system is not altered as a result of adiabatic 
changes. In agreement with this postulate, the quantum 


conditions would be stated as 
nh = f(z, Q1, Ae, ae) ) 
where f is an adiabatic invariant. 


a. Properties of adiabatic invariants 


Let the Hamiltonian function be H(p, q; a, d, .. 


the a; are adiabatically variable parameters: 


.), where 


(The a; do not enter here because they are supposed to 


be small.) Therefore, 


aH ab _y (2H) , 
ap dt Ver. 


For periodic motions 


Since f is an adiabatic invariant, 


df/dé = 
and 
df of dH 
dt OF dt a Sha, 


therefore, the conditions for adiabatic invariance are 


(se) (Se), (Sa) = 
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b. Examples: 


1. Linear oscillator 


pm 
Tomes oH _ ; DB oot 
Ow @ , 
a=, ag =m 
= oH oH = = 
Seed ae laa eet Ein = Eno = $B 


Therefore, 


It is easy to verify that f=2#/w is an adiabatic invariant: 


on? 
om ~—_ 2m? 


1 
= gore = a (— Exin+ Epot) =0. 


2. Nonrelativistic particle on a line 


H=p?/2m. 
The force on the wall (Fig. 21.1) is 


sy 2S 
Therefore, 
Hl? = constant. 
l 
fixed 


movable 
Figure 21,1 
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OZZZLLLL LLL 
— (p+dp) 
Figure 21.2 
Varying the mass, we obtain 
SH = — ESm/m 

or 

Em = constant . 
Therefore, 


Emil? = constant . 


The energy shell of the system is given by (see Fig. 21.2) 
dQ=2ldp. 


However, for nh=2l|p|, the density of states in phase 
space is 1/h. Therefore, 


From the previous considerations, this is indeed adiabati- 
cally invariant. The most general way of satisfying both 
conditions is obtained by replacing n by f(n), where 


eo ae 


uc nN 
3. Relativistic particle on a line 


K=2pv/2l, 


vdp = 8E=—K31=— 2p dl, et Aeon 
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Therefore, 
pl= constant . 


The conditions for adiabatic invariance and the condition 
on the density of states are still satisfied, because 2l|p|=nh. 
Only the energy eigenvalues change: 


22 
Ib). = Vm pn : 


In the limiting case m—0, we have v=c and E/p=c. 


Therefore, 


nh 
Ei, = ¢ > 


(as for the radiation cavity), and 
Hl = constant. 


For radiation, 


By comparison with the limiting relativistic case, 


h 
p= and B= =h. 


These equations were the point of departure for de Broglie. 


22. VAPOR PRESSURE. NERNST’S THEOREM 
For a monatomic gas, we have (see Section 20) 
w= kT {log p— $ log T—7}4 Ky, 


7 2am)? k? 
png CR 


where H, is the zero-point energy per atom. 
For a diatomic gas, we can no longer calculate classically, 
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because the rotational degree of freedom about the sym- 
metry axis must usually be considered to be frozen in. On 
the other hand, at room temperature we may usually cal- 
culate with the other two degrees of freedom classically. 
In that case, we have per molecule: 


1 teal 
Vie az [exe |- kT 2A (xi + 7)| dz, dz, sin§ d0 dg 


For tri- and poly-atomic gases, we obtain 


al if 
Za 55 [ex| — an +34 72) aan dnd, sin 9 a0 dpdy 


1 
= 5, (2nkT)#( A, Ay Aa)t X dex 2m. 


For molecules with two or more identical atoms there is 
another complication. To illustrate this, we consider dis- 
sociation equilibrium in a gas. Let the molecule be of the 
type M=(C,),{C,), ; that is, let it be composed of p, atoms 
of type OC, and p, atoms of type C,. Let there be N mole- 
cules of type M, N, free atoms of type C,, and N, free atoms 
of type C,. Let C, be the total number of atoms of type C,, 
and C, be the total number of atoms of type C,; C, and OQ, 
are fixed numbers: 


Ni+pN=C,, 

Nz + piN = C,, 
exp[— F/kT] =exp{— Fy/kT] exp(— F,/kT] exp[—F,,/kT]P . 
The factor 1/h’ is to be thought of as already incorporated 


into exp[—F/kT]. The combinatorial factor P is equal to 
the number of different systems that can be formed hold- 
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ing N,, N,, and N fixed: 


p= number of permutations of identical atoms 
number of internal permutations 
SECs 
NIN! No" 


The o which appears in the denominator is called the sym- 
metry number. It is the number of permutations of identical 
atoms which can be generated by rotating a molecule. In 
other words, it is the number of equivalent orientations of 
a molecule or the order of the symmetry group of M. For 
example, 


Gi. o— 1, 
H,: «o=2, 


Letting i=1 and 2, and defining 


1 
C Wap Oto cxPl- F/kT}, 


exp[— #*/kT] = 
i 
exp[— #4/kT] = wicse [— F,/kT], 
and 


exp[— Fu /kT)= =, exp [— FalkT], 


at 
N!o% 
we obtain 
exp [— F*/kT] = exp[— Fy,/kT exp [— Fi/kT]exp[—F? [kT]. 
Applying this to our problem of the vapor pressure, we 
obtain the same yu for the polyatomic gases as for a mon- 
atomic gas, except for the following two modifications: 


SlogT— (C,/R)logT, where 


= (monatomic), 
les: =+% (diatomic), 
e =4 (triatomic). 
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2. The quantity i is different: 


w=kT {logp—(1+ $f) -logT—1} +E, 


; (220m)? k? 
monatomic: = -log oC we 
; ; (2am)tk? 82 Ak 
diatomic: i=log — ao 
2am)*k* 87? (27k)? (A,A,A;)? 
tri- and poly-atomic: i=1lo (| ” mle es 2 ) 


a. Applications 
1. Solid-gas vapor pressure curve: Using 


v 


: 9 hy 
Peous = Ey+ kT = 5 |e dy log (1 — exp | im) 
o 


along with the initial values, we can calculate the curve. 
2. Chemical equilibrium in a gas: 


W*(N,, M2, W) 
= constant x exp |- a [F*(N) + FI(N,) + P20] , 


For the most probable distribution, we obtain 
SlogW*=0 for 35VN,=—p,8N and 8N,=—p,SN. 
Thus, 


1 8Ny+ oS Na+ Wy SN = 0 
or 


Lu Pifi— Pel, = 9, 


which is the result obtained with thermodynamics. As a 
consequence, we now know the constants that appear in yw. 
Furthermore, it is also possible to calculate fluctuations in 
this manner. 
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b. Quantum theory and Nernst’s theorem 


Nernst’s theorem, which says that all bodies at absolute 
zero have the same entropy (normalizable to zero), implies 
something about the density of quantum states in phase 
space. We have 


F=—kTlogZ, and Z=)4g,exp[—£#,/kT], 


where a degeneracy g, for state n has been assumed. If 
the system is finite, then there are only discrete energy 
eigenvalues. Let E, be the lowest eigenvalue and EF, be 
the second lowest. For k7?<EH,—E,, we can make the fol- 
lowing. calculation: 


.. 7 I _ £,—£, 
7 = exp|: Blk T}(1+ exp] iT J+), 


ms E, 9; E,—E£, 
log Z = — 74 + loggs + Lexp |— iT |+-, 
E,—E 
F = By— kT log go— kT exp |- a l, 
oF E,— E,\ g E,—E£. 
S= ——=kl k pa oD) Ve 2 0 ; 
Or ogg + ( ae ) Sexp | kT 
Letting T—0, we obtain 
S=hklogq. 


Therefore, Nernst’s theorem requires g,=1. This means 


1. The lowest state is not degenerate; 
2. The lowest state is sufficiently separated from the next- 
to-lowest state. 


23. EINSTEIN’S THEORY OF BLACKBODY RADIATION 
Shortcomings of the method of oscillators are 


1. E,= (c?/82r*)o, is a classical formula; 
2. Molecules are not oscillators. 
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We now assume an arbitrary system with energy eigen- 
values E,, H,,..., but we do not specify either selection 
rules or degeneracies. On the other hand, we use the Bohr 
frequency condition: 

LY, — B= Wee. 


Einstein > introduced transition probabilities which describe 
the behavior of the system statistically: 


spontaneous emission n—>m: Aj, 
(E,>E,,) }induced emission > > Ba Os, 
absorption m—>n: By oy. 


The number of transitions in a time interval dt is given by 
n—>m (At, -+ Beopaladt = dZeg5,,., 
m—>n Bo 0y Nin dt = AZ mon 
Now assume thermal equilibrium: 
N, = 9.0 exp[— E,/kT}, 


AZ nom = AZmon « 
Therefore, 


gn(Am + Brgy) exp[—E,/kT] = 9, Bro, exp[— E,/kT], 
Apa sos ( Bee ce exp[(Z,— B.)/kT)) : 


Substituting (Z,—Z£,,)/h=y», we obtain 
= A;,/B*, 
oF Gn Be [gn Ba) exp [hv] kT] — 1" 


If, for high temperatures, we require that this reduce to 
the Rayleigh-Jeans law (see p. 71), then 


(1) IuDa — Gaba ) 
(2) = a sl 
an C 


5 Physik. Zeitschr. 18, 121 (1917). 
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For arbitrary temperatures, Planck’s law follows from this 
immediately. The problem of calculating the A”, B", and 
BY separately can only be solved with quantum mechanics. 


Momentum transfer between radiation and matter 


We consider the Brownian motion that a particle under- 
goes in the radiation field. From the results of Chapter 3, 
we have 


d.=|Xat, G,=0, G2=2kTWr. 
Oo 


1. Calculation of the resistive force W acting on a particle 
moving slowly in the radiation field. Denote the rest system 
of the radiation by K and the rest system of the particle 
by K’. During absorption, the recoil momentum is, for 
the transition n—m in K, 


where 9 is the angle between the photon and particle mo- 
menta; during emission, the recoil momentum is 


ah cosé. 
C 
Define 8 by 
S = 9, exp[— #,/kT] + gn exp[— E,,/kT]+-... . 
The fraction of time which the particle spent in state n is 
Jn ©XP [—E,,/kT\/S ; 
the fraction of time spent in state m is 
In exp(— E,,/kT]/S8 . 


The numbers of absorptions and emissions into solid angle 
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dQ’ are, respectively, 


dQ’ 
59m exP[— En/kT] B06 


and 


il q d.f2" 
‘ Bre : 
5 9n exp [— Bp /kT}Br 06 


(Note that the prime refers to the system K’ and that o, 
is not isotropic.) Let » = (H,—E,,)/h. Then, 


ye 


= gn Be (exp[—Em/kT]— exp[— B,/kT) 


t 


ss ,dQ 
x {6 ) cos 0 ys 


Without proof, we state that the following relation holds 
[A-7]: 
Qr _ & 
y's ys ° 
Further, to first order in v/c, we have a Doppler shift given by 
v 
y=» (1 a 7 e036) - 


Therefore, 


' 0 . 
ey = (. ae - cos a’ O[1+(v/0) cos O' Jy’ 
a (1-32 coss')] st (= ~ cos 6’ 
Ov, = C Oy, a se iF cos | 
~ oy +0086! Og 3 
— Ov, C Vo ny = Ov, | > 
; dQ’ v 1 (do, 
fame ——tfa in) 
er 
—  ¢@ Ov \v3] fy,” 


Then, again writing » instead of », and interchanging n 
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and m, we obtain 


hy 1 


=— yo By exp[— E,,/kT](1 — exp[— hv/kT7]) 


[350] 


In equilibrium, we have the Planck formula, 


82th ys ; 
Oo? exp[hv/kT]—1’ 
7 (*) _ 8ah — exp[h»/kT] h 
e} (exp[hy/kT] —1)? kT 
~.! h exp [hv/kT] 
1 ET explhv/kT]—1 


1 fhe? 1 g 
VY=-~(—| —= — E,'k sg 
V 5 (2) Sel ERA BS 0- , 


(Gm/S) exp(—#,/kT| Bo, = number of absorptions per se- 
cond = 4Z (Z= total number of all processes per second). 


Hence 
1 fhve\? 1 
eee (7) onT 2 


2. Calculation of G*. 


where Zr= total number of all transitions (induced and 
spontaneous emissions plus absorptions) during a time t. 
Interpretation. Each transition, including spontaneous 
emission, imparts a recoil momentum hy/c to the atom: 
G? = (2) cosas = (2): 
Cc 3\e 


Thus, we are compelled to ascribe a recoil action even 
to the spontaneous emission process. This is in contradic- 
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tion to interference observations, from which it is known 
that spherical waves can be coherent. 

This paradox is resolved by the remark that an exact 
knowledge of the momentum of the particle is required for 
a@ measurement of the recoil momentum. According to the 
Heisenberg uncertainty principle,* exact knowledge of the 
Momentum would exclude the exact knowledge of position 
that would be necessary for observation of interference. 


24. QUANTUM STATISTICS OF IDENTICAL PARTICLES 


In the case of radiation, the state was completely deter- 
mined by giving the energy n,hy, of each normal mode. 
In particle language, this means giving the number of 
photons n, in state s. One can try to use the same de- 
scription in the case of an ideal gas. Let the word ‘“‘state”’ 
denote the state of an atom and the word ‘‘cell’’ denote 
the state of the entire system. According to our previous 
considerations (see Sec. 7), for a specified system (n,), there 
are 

n! 
= M1!Mg!... My... 
distinct cells. ~ 

If we consider radiation, then, because of the wave nature 
of radiation, the situation is different. Guided by the wave 
nature of matter, we can tentatively try the same procedure 
with gases that we used for radiation. 

For low densities, such that n,—0 or 1 for all s, we 
have P-—-N!. The density of states is then 


1 


win® instead of 


1 

* . 

In the particle picture, this has the consequence that the 
particles are no longer statistically independent. For ex- 
ample, consider a case for which there are two states, A 
and B, and two particles. 


* See W. PAuLI, Wave Mechanics. 
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Statistical Weight 


ns Np L.E. statistical independence F.D. 
2 1 i 0 

1 1 1 2 ] 
0 1 1 0 


In the case considered here, the ‘‘statistics of symmetrical 
states’ or Bose-Einstcin statistics (B.E.), the particles have 
the tendency to condense into groups. 

In addition, there is the ‘‘statistics of antisymmetric states’’ 
or Fermi-Dirac statistics (F.D.), for which the n, are res- 
tricted to have only the values 0 and 1. In this case, it 
is evident that the particles repel one another. 

In the case of radiation, the total number of particles is 
not specified. In order to preserve the analogy with radi- 
ation as nearly as possible, we start from the grand canon- 
ical ensemble (see Sec. 14). Then we do not have to write 
the N! because the particles are now considered indistin- 
guishable. 


(nr, nay... oe, «..) = exp [alO + eN — Bimy, ...5,...)}] - 


For radiation, the w in this formula must be set equal to 
zero, since WN is not an independent variable. 
We now consider an ideal gas. Let ¢,=the energy in 
cell s. Then, 
D> Nee., 


v= i, 
and 


W*= exp [a2] [[ exp[a(u— e,)n,| : 


From > W*=1 we obtain 
all 


— exp [ou = é,)]’ 
De 1, — Ori: 1 = exp[aQ] TI{1+exp [(u—e,)]}. 


met 2, = 0, ...,. 00% 1= exp[a2] TT; 
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Therefore, 
B.E.: aQ= > log {1 — exp[a(u — e.)]}, 


B.D.: a0 — Soe (exp aly —,)]}. 


From N=—02/du we obtain 


i 
7 explo — 2) — 1 
al 


BD N= 2 expla el) Fl 


Bew.: = >} ———_ 


From the requirement N>0, we see 


B.E.: —comcapn< 
ED ee 


From §=—(02/0T), we obtain 


B.E.: a eel exp[a(u— €,)] 

a. | a sa 
B.E. gu —F+E —2-ypN+E 

F.D. a rT, ? 

where 


Es 
i 2 pie e(a— el 


a. Transition to integrals 


For an ideal gas of mass points the density of states is 


dZ _y 27 (2m)? 


de h3 a 


=o =a, de=kTdx, and A = op = 
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we obtain’ 


Q= -vEn spbiomit 


ner F exp[A— 2]})V/@de. 
0 
With p=—(0Q/OV),,, we obtain 2=—pV, as it should 
be. Integrating by parts, we obtain 
+ flog (1 = exp[A — 2}) $d(c) 
Qo oa 
= x §at log(1 = exp[A — 2])}— 


5 fe exp[A— a#]dx 
3 1—-exp[A—a]’ 


and hence? 


02>, vag? i pee 


he 3Jjexp[—A+a]F1° 
0 
Further, 
Qn ame r Bae 
Hees kT) | ——_—_—____ 
eee J exe A+«]¥1’ 
22 (2m)* j atdxr 
tal kT) : 
h3 — exp[—A+2a]+1 
Therefore, 


—Q= pV=2E/3, 


which is the result for an ideal gas. 
With the definitions 


2 { VEdx 
Va J exp[—A+2]F1 
5 


F;(A)= 


and 


@ 


atdr 
— a ne A+a|]=1 


7 For Bose-Einstein statistics, the upper sign is to be taken; for Fermi- 
Dirac statistics, the lower sign is to be taken, 
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we obtain 
2QmnkT)* 
N= y! ae F(A), 
3 (2am)! (kT)? 
= Vis —— Gz(A). 


Here, A=yp/kT is to be thought of as a parameter which 
is to be eliminated from the two equations. By partial 
integration we find 


b. Limiting cases 

1. Dilute F.D. or B.E. gas (A——oco). 
if exp[—|A|—@] 

exp[—A+a]+1 1-Fexp[—|A|—<2z] 


= $(41) exp[(—|4|—a)n], (4<0). 


n=] 


With 


@ 


aa 7a joel na \etda = Ea 


and 


t-<) 


alt 
= = [exp [—nax)xida = re 
0 


we obtain, for A< 0, in the B.E. case, 


FocAn eellet . 
C= pee 

and in the F.D. case 
F(A) = x yy PL = bay 
G4(A) =F (-1y SPL A 
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For |A|>1, only the first term in the sums need be con- 


sidered: 
2mnkT)* 
yoy sues = Ve, 


Sau V (2makT)t 
A= = log | Ome : 


This is the w of a normal monatomic gas. Further, from 
C2 we 


follows 
E=8NkT, —2=pV=NkT. 


To a higher approximation we obtain per mole 
L 
oV=RhL (: Slr Ns: Tv (Q2amkT)-* + is) : 


2. F.D. gas at low temperature (A—oo). For F,(A) and 
G,(A) there are the following asymptotic series: 


Bild) = poe Al (14 +.) 
5 2 
GA) Tee (1 +epte). 

Proof: 
Let 

_[__fa)an =A 

-|— “tee % TY; 
0 


f(A+y)=gly), and Jor ay =a. 


—A 


Then, 


G(y) evdy 
(ev+ ily Ps 
—A 


_fgydy _ ( eyydy _ 
Fett J ew+1 e+ 
0 


| —A 
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The Taylor series for G (with the remainder) is 
G(y) = G(0) + y@'(0) + 3y7G"(0) +... . 


Hence 


+0 


i a ST) | ee, dy 


—o 


+a 
; POG aay yer dy 
sas fee TR cs 


—o —o 


T= G0) +0) +..., 


@(0) = [aor a0= fr de, @'(0) = f(A), 60) =f'(A), 


f(x) da r a 
=f alt ies =|reyae+ = (A+ 


From this the assertion follows immediately. 


Now, let f(a) = f(e/kT) =9(e), and let A=y/kT. Then we 
have 


td 
(e) de or z fy Ad -, 
ler maa =} ple) de + g (ft) p(t) +.... 


In the limit as T-—>0, the integral on the left can be re- 
placed by the first term on the right, which is an integral 
whose range only extends to w. The function is given 
by p(e)~et for N, and ~e! for E. 

Zeroth-order approximation: 


27 (2m)! 2 a 2 
h3 3 Ho , =v 5 ua : 


N= Vv 


The zero-point energy becomes 


m= ¥ als) mr) 
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First-order approximation: 


wey 2 eae ary| 


hs 


3 
27 (2m) 2 
5 


i he 


5702 
ap! ar (| 
pla ) 


3 70? 
= Nop \1+ — ahi, 
=e | +5 (er 
If we keep the definition of ~, used in the zeroth-order 


approximation, then we obtain 


702 
ul E + 38 ary| = mg 


or 
4 


2 
i, ji aa er +-..| 


3 5a 
H=N— 
Nghe i+ 


(kT)? + | ; 


2 
0 


OH 70? 
aes | ) cea pen /-0 
C, (ar) 5a" ee 


For the entropy we have 


_E-F _E—yN+pV _$E—~N 


2 uly i i i 


2 
Sy ap (Nernst’s theorem) . 
Zio 


Application to an electron gas: In order to take account 
of the electron spin, we can consider the gas as a mixture 
of a two-component gas with opposite spin directions. In 
the absence of external forces, each gas contains N/2 par- 
ticles. Therefore, the previous formulas are valid if we 
replace N by N/2 and EL by H#/2. For u.>kT, we obtain 


Ey — Niu , 
27 (2m)* 


N=2V = 


Bit 
3 


O° 
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In the presence of an external magnetic field, there is 
an additional energy u,H, where u,= eh/(2mc) is the Bohr 
magneton. Let us denote the cases of spin parallel to H 
and antiparallel to H by the subscripts 1 and 2, respec- 
tively. The equilibrium condition is 


by fig = M2 + fel 
This equation and the equation 


27 (2m)? 


Nv) i 


(Zui + 3ui) 


together determine yu, and w,. Further, we have 


M= pa( N, — N,) = baV 


and 


This is a weak, temperature-independent paramagnetism. 
(The result is only qualitative, since forces between the 
electrons have been neglected.) Superimposed on this pa- 
ramagnetism is an orbital diamagnetism which is } as 
strong [A-8]. 


3. B.E. gas. Degeneracy. The series given for the first 
limiting case are always convergent for the B.E. gas. 
Particularly, for d=0, the values are F(0)—2.615 and 
G(0) =1.34. Hence 


V 
vs, = hs (QamkT)i °2.615 = Nee ? 


1.34 
2.615” 


1.34 
2.615 © 


Ey = 3N kT 


= $hN, 
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That is, for fixed V and 7, there is a maximum number of 
particles that can be in the gas. If we try to add more, 
then, following Einstein, what happens is the following: 
The particles go into the condensed state («=0, S=0, 
p=0, E=0, and, therefore, w=0) and, in a sense, form 
a second phase which is in equilibrium with the first phase, 
since w= 0 for both. 

One can object that replacement of the sum by an in- 
tegral is no longer correct for A=0. If the calculation is 
carried out with sums, then this bounded WN, is not ob- 
tained since for ~=0 the function is singular at e=0: 


JE ma i 
exp[e/bT]—1 = \/e) 
The result depends on the type of vessel. If the vessel is 


sufficiently large, then the above theory is nevertheless 
essentially correct. 
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Appendix. Comments by the Editor 


[A-1] (pp. 12, 14, 15, 17). The Boltzmann equation is usu- 
ally written in the form 


a (ai). 


df zs Of. 
file Ox 


where 


of 


is the kinematic term and 


(a) =| favanyo WV) — f(w)f(V) wg 


is the collision term. There exist two distinct situations of 
importance where the collision term is negligible. 


1. Local equilibrium: (Of/Ot)...= 0. By definition a local 
equilibrium distribution annihilates the collision term. This 
implies that Eq. [4.4] on p. 10 must hold. Now this equa- 
tion as well as Eqs. [4.4a], [4.7a], [4-7b] are not changed 
by letting f depend on x and ¢ and by including a poten- 
tial energy #,,,. Therefore the most general form of a 
local equilibrium distribution is 


flv), i) = Aexp{— B(x, o|F 38 + Big, t) — mu elas all. 


where any (x, t)-dependence of A can be absorbed into £,,, 
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which therefore may also depend on time. The function 
B-+(x,t) = kT (x,t), where k is Boltzmann’s constant, de- 
fines a local temperature T(x, t), and c(x, t) is a local drift. 
The distribution f, is restricted by the condition that the 
collision term describes the collisions properly which means 
that the frequencies involved are small compared to the 
collision frequencies. In addition, it has to satisfy the 
Boltzmann equation, which now reduces to the kinematic 
equation 


A stationary distribution is a special case of an f, because 
of/ot = 0 implies d#/dt = 0 and hence again (0f/0t),., = 0. 
However, now the Boltzmann equation reads 


of of 


an Ut ay Kim=0, 


which implies that E,,, is time-independent and B= 2a/m = 
= constant; c = c°+w x xis perpendicular to K = — 0F,,,/0x 
and c° = constant, w = constant. 


2. High-frequency limit: |Of/dt|>> |(Of/Ot)on|. The Boltz- 
mann equation again reduces to df/dt= 0. A high-frequency 
solution of this equation is called a collisionless distribution. 
This case is of importance for high-density Fermi liquids, 
in particular, electron plasmas. However, because of the 
strong interparticle forces, one has to include in the po- 
tential energy a distribution-dependent term of the form 
fd?’ p(v, v')f(v’, x,t), which gives rise to a so-called Vlasov 
term in the collisionless Boltzmann equation and to a col- 
lective mode (plasmon, zero sound). 

An approximation of the collision term often used for 
transport processes (see Sec. 5) is the collision time appro- 


imation, 
= 
(f) =-U-tor, 
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where f, = A exp [— a(v — c)*] is a stationary distribution. 
Ii f— fy oscillates as exp (— iwt) and if the system is per- 
turbed by a weak external force K = K,,,, the solution of 
the Boltzmann equation becomes, to lowest order in K,,,, 


T 2a 
f= h+ia WK} ho, 


where u=v—c. Note that a temperature gradient may be 
considered as an external force according to the substitution 


2a m 1 cee 
== — 
Tee” yer Oe 
with 2«/m=1/(kT). For w= 0 the above f then becomes 
a special case of the expression given on p. 17. 
The frequency range wt <1 is called the hydrodynamic 
domain; wt>1 is the collisionless domain. 


[A-2] (p. 29). With modern computers it has become 
possible to calculate such time averages and hence, in prin- 
ciple, to check the ergodic hypothesis ‘‘by experiment.”’ 
See, for example, A. Rahman, Phys. Rev. 136, A 405 (1964); 
L.:Verlet, Phys. Rev. 159, 98 (1967). 


[A-3] (pp. 32, 39, 53, 54). The virial theorem 
2NEy, — (V™"+V™) = N ——-q?=0 


follows trivially from Newton’s equation and the formula 
for the Brownian motion q*=2Dt derived in Secs. 16, 
17, 18. Here 


| fa Sa be = pq-dS = 3pV 


is the external virial due to the forces acting on the walls. 
For central two-body forces 


F** = _ > qi ~ ONG 3)5 yy = lq: a qi| ’ 


eee Vik 
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and the internal virial is 
int 
yixt a Sak. a 4 Vik U' (rx) 


By partial integration one finds per mole, assuming that 
U(r) tends to zero stronger than r-* for 7 -> 0, 


© 


ss 2 
AT) = = L?o3 U(a) — 3 fe 4or*dr . 


Now U(c)=0 by definition (see Fig. 13.3), but application 
of the virial and equipartition theorems (p. 32) to the radial 
motion gives 


U(c) = —rU (r) = 2B = kT. 
Together with Ey = 3kT the expressions of p. 54 then 
readily follow. 


[A-4] (p. 56). The functions Ff, F,, and F, are the same 
because homogeneity implies that the Hamiltonians are 
sums over the respective particle numbers N, N,, and N, 
of the same 1-, 2-, etc. particle Hamiltonians. 


[A-5] (p. 57). This ceases to be the case in the vicinity 
of a phase transition where the kth Fourier component of 
the fluctuations, (AN,,)?, is of the form 


a constant 
- (AN,,)? = (E24 kya Te k)i-n/2 * 
ere 
eT) = constant 


. (‘a a 


is the correlation length and T, the transition temperature. 
While classical theory (Ornstein-Zernike) yields for the crit- 
ical exponents vy = 0.5, 7 = 0, experimental evidence favors 
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slightly larger values of these parameters. For more detail 
see, for example, H. E. Stanley, Phase Transitions and 
Critical Phenomena (Oxford University Press, New York, 
Loi). 


[A-6] (p. 85). The meaning of the zero-point energy which 
Planck had introduced somewhat reluctantly in 1911 was 
a much discussed. problem in the days of the old de Broglie— 
Sommerfeld quantum theory. Pauli discussed it extensively 
with Otto Stern in Hamburg in the early 1920s. Stern had 
calculated, but never published, the vapor pressure differ- 
ence between the isotopes 20 and 22 of Neon. Without 
zero-point energy this difference would be large enough 
for easy separation of the isotopes, which is not the case 
in reality. 

Pauli, on the other hand, had calculated, but also not 
published, the zero-point energy of radiation and came to 
the conclusion that its gravitational effect would be so 
large that the radius of the universe ‘‘could not even reach 
to the moon”’ (short wavelengths were cut off at the clas- 
sical electron radius). 

For more details on this question see C. P. Enz and 
A. Thellung, Helv. Phys. Acta 33, 839 (1960). 


[A-7] (pp. 85, 96). Let the adiabatic compression of the 
cavity of p. 85 be produced by a piston moving with ve- 
locity v. Then a photon of frequency » incident on the 
piston under an angle 6 will be reflected with frequency 1’ 
determined by the Doppler shift formula 


v= v(1-+ = 0088) 


The work done by the radiation pressure p = J,(2/c) cos 6, 
where J, is the intensity of the radiation, during a time 
interval St of compression is per unit surface, 


pvdst = 12 cos 68t = (1, — I,)8t . 
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It follows that 

La (1 + cos ) ; 
and hence . 


IR 


dr 
v 


But EF, is proportional to Z,, which proves the relation 
quoted on p. 85. 

For the particle on p. 96 the Doppler shift formula does 
not contain the factor 2 because of absorption instead of 
reflection. But otherwise the situation is the same since 
the particle is supposed to move slowly in the system K, 
so that the transition to K’ is adiabatic, that is, does not 
induce Bremsstrahl photons. Therefore the relation £,/y 
= Ej, /v' of p. 85 also holds here. With the expression 


wee 
» gee CP 


of Sec. 19 (p. 71) the relation quoted on p. 96 readily follows. 


[A-8] (p. 106). The formula derived for WM is called the 
Pauli paramagnetism; see W. Pauli, Z. Physik 41, 81 (1926). 
The orbital diamagnetism which is due to the Larmor pre- 
cession of the electrons and has the value — 4M is called 
the Landau diamagnetism; see L. D. Landau, Z. Physik 64, 
629 (1930). 

In a real solid the electrons move in energy bands and 
the formula is modified. The diamagnetism of single-band 
electrons was calculated by R. E. Peierls, Z. Physik 80, 
763 (1933). The multiband case is considerably more com- 
plicated. For a review see C. P. Enz, in Semiconductors, 
Varenna Summer School 1961, edited by R. A. Smith 
(Academic Press, New York and London, 1963), p. 458. 
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